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PREFACE 


This  report  serves  as  a  user’s  guide  to  present  the  start  up  procedures  and  general  formulation  of  the 
Articulated  Total  Body  (ATB)  Model  Version  V,  a  dynamics  simulation  tool  for  the  human  body 
biomechanics  in  various  dynamic  environments. 

To  provide  a  complete  user’s  guide,  Obergefell,  Gardner,  Kaleps,  and  Fleck’s  “Articulated  Total 
Body  Model  Enhancements,  Volume  2:  User’s  Guide,”  AAMRL-TR-88-043  (1),  has  been  modified 
and  incorporated  into  this  report. 
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1.  INTRODUCTION 


The  Articulated  Total  Body  (ATB)  Model  is  used  by  the  Air  Force  Research  Laboratoiy  (AFRL)  and 
other  organizations,  companies,  and  educational  institutions  for  predicting  gross  human  body 
response  in  various  dynamic  environments,  especially  automobile  crash  and  aircraft  ejection  with 
windblast  exposure.  The  ATB  Model  originated  from  the  Crash  Victim  Simulation  (CVS)  program 
(2).  Aerodynamic  force  application  and  a  harness  belt  capability  were  added  to  the  CVS  program 
by  Calspan  Corporation  in  1975  for  AFRL  (3),  and  the  resulting  program  became  known  as  the  ATB 
Model.  In  1980,  Calspan  made  a  number  of  modifications  to  the  ATB  Model,  combining  it  with  the 
then-current  3-D  Crash  Victim  Simulation  program  to  form  the  ATB-II  Model  (4).  Complete 
documentation  of  the  program  through  the  ATB-II  version  was  performed  by  Calspan  (5,  6, 7, 8). 
The  next  version,  ATB-in,  which  included  improvements  made  by  J&J  Technologies,  Inc.,  was 
generated  to  model  the  body  response  to  windblast  for  AFRL  (9).  The  version  ATB-IV  (1, 10, 1 1) 
was  released  in  1988  with  a  number  of  additional  efforts  being  made  to  improve  various  aspects  of 
the  ATB-m  Model,  with  emphasis  on  its  capability  to  simulate  aircraft  ejection  with  windblast 
exposure,  as  well  as  complex  automobile  accidents. 

The  ATB-V  Model  introduces  three  new  simulation  tools:  water  force  simulation  (12),  joint 
actuators  (13),  and  deformable  segments  (14).  A  major  change  has  been  made  to  the  data  arrays  to 
increase  the  maximum  number  of  segments,  planes,  and  contact  definitions.  A  new  type  of 
structured  ASCII  graphics  data  output  file  has  been  designed  for  use  by  the  Interactively  Manipulated 
ATB  Graphical  Environment  (IMAGE)  program  (15)  and  the  VIEW  program  (16).  It  has  a  more 
efficient  format  that  the  original  graphics  data  file  and  mikes  troubleshooting  easier  for  the  IMAGE 
and  VIEW  programs.  Several  outdated  features,  such  as  the  restart  and  plotting  options,  have  been 
deleted  from  the  ATB  Model.  There  are  also  a  number  of  other  minor  modifications  and 
input/output  enhancements,  such  as  enhanced  HIC  value  computations,  clarified  joint  property 
definitions,  calculation  of  moments  of  inertia  and  principal  axes  for  a  set  of  segments,  and  control 
of  individual  contact  ou^ut. 

This  User’s  Guide  accompanies  the  release  of  the  ATB-V  version.  It  contains  comprehensive 
information  on  the  ATB  Model  and  its  input  structure.  It  is  completely  restructured  from  the 
previous  version’s  User’s  Guide  with  extensive  modifications.  Section  2  gives  a  general  description 
of  the  ATB  Model  and  its  structure.  An  overview  of  the  ATB  input  data  and  ou^ut  files  is  provided 
in  Section  3.  The  appendices  contain  example  input  and  output  files  from  the  model. 

1.1  Installation  and  Hardware  Requirements 

The  ATB  Model  is  written  in  FORTRAN  so  it  can  be  run  on  a  number  of  platforms,  including 
personal  computers  and  UNIX  workstations.  The  software  package  includes  the  FORTRAN  source 
code,  the  PC  executable,  and  example  simulation  files.  For  installation  information,  please  see  the 
file  readme. 
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The  PC  executable  requires  a  486  or  higher  personal  computer  (PC)  supporting  MS-DOS  mode 
operation.  It  is  recommended  that  the  PC  have  at  least  8  MB  RAM  and  5-10  MB  free  hard  disk  space. 
The  program  is  typically  distributed  on  3.5-inch  high  density  (1.44  MB)  diskettes. 

1.2  Running  the  ATB  Model  —  An  Example 

In  order  to  demonstrate  the  procedure  for  running  the  program,  the  example  file  shipped  with  this 
program  is  used  as  the  input  file  in  the  following  discussion.  It  is  listed  as  cxctmplcMin  on  the 
shipment  disk  and  is  presented  in  the  Appendix.  This  example  simulates  a  sled  test  with  a  male 
subject  in  a  seated  position  facing  forward. 

1.2.1  Programs  and  I/O  Files  in  the  ATB  Model 

The  executable  program  of  the  ATB  Model  is  named  atbv-x.exe.  There  are  two  other  executable  files 
shipped  together  with  the  ATB  Model.  One  is  the  GEBOD  program,  gebodx.exe  (17),  which  is  used 
to  generate  the  human  and  dummy  data  sets  for  the  ATB  input  file,  i.e.,  the  B  cards  of  the  input  file. 
The  other  is  the  VIEW  program  viewx.exe  (16),  which  is  a  graphics  program  used  to  generate  the 
pictures  of  body  motion  based  on  the  simulation  results.  An  alternative  graphics  program,  called 
IMAGE  (Interactively  Manipulated  ATB  Graphical  Environment),  is  available  on  Silicon  Graphics 
workstations  to  generate  solid  object  animations.  Figure  1  presents  the  organization  among  these 
programs  and  their  corresponding  input  and  output  files. 

When  using  the  ATB  Model  for  occupant  modeling,  one  typically  uses  the  GEBOD  program  to 
generate  the  occupant  segment  and  joint  data.  The  GEBOD  program  creates  the  B  cards  for  ATB’s 
input  file  and  saves  them  in  a  file  with  an  .ain  extension.  The  generated  B  cards  are  inserted  into  the 
input  file,  which  also  has  the  extension  .ain.  (The  ATB  Model’s  input  file  requires  an  .ain  extension 
following  the  file  name.)  In  the  example  case,  it  has  the  full  filename  example.ain.  Upon  successful 
execution  of  the  program,  a  set  of  output  files  is  generated  with  a  user-chosen  filename  and  a  set  of 
unique  extensions  to  distinguish  them.  There  are  three  kinds  of  extensions  assigned  to  the  output 
files,  .aou,  .tpl  or  .sal,  and  .t??.  The  question  mark  stands  for  a  number  character.  The  file  with  the 
.aou  extension  is  the  main  output  file  containing  all  the  input  data,  clearly  labeled,  and  extensive  run 
time  information.  The  file  with  the  extension  .sal  or  .tpl  is  the  data  file  used  by  the  VIEW  and 
IMAGE  programs  for  graphical  image  generation.  The  set  of  files  with  .t??  are  the  tabular  time 
history  outputs  of  user-chosen  simulation  data,  such  as  segment  accelerations,  contact  forces,  joint 
forces  and  torques,  etc.  The  number  of  time  history  files  is  defined  by  the  user  input. 
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Physical  occupant  data 
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Control  file  ^ 
I  *.setup  file 

Figure  1.  Organization  of  ATB  Program  and  Files 


1.2.2  Running  ATB  on  PC 

To  start  running  the  program,  make  sure  you  are  in  MS-DOS  shell  or  mode  and  type; 

atbv-1 

then  press  [Enter].  You  will  be  prompted  by  the  message  “ENTER  INPUT  FILENAME, 

EXTENSION  .ain  IS  ASSUMED:”.  To  run  the  input  file  example.ain,  type  in: 

example 

then  press  [Enter].  You  will  be  prompted  by  another  message  “ENTER  FILENAME  FOR  ALL 

OUTPUT  FILES,  EXTENSIONS  WILL  BE  ASSIGNED:”.  Type  in: 

example 

then  press  the  [Enter]  key.  You  can  choose  another  output  file  name  instead  of  example  if  you  wish. 
The  program  will  start  ranning  and  generate  the  output  files.  If  no  error  occurs,  the  program  will 
terminate  automatically  with  a  message  “COMPLETED  ATB  SIMULATION”  displayed  on  the 
screen. 

The  ATB  Model  also  has  a  workstation  version  ranning  under  UNIX.  The  ranning  procedures  are 
exactly  the  same  as  for  the  PC  version. 

1.2.3  Runtime  Error  Message  and  Troubleshooting 

There  are  two  types  of  runtime  error  messages.  One  is  from  the  compiler  or  operating  system.  For 
these  errors,  consult  the  relevant  reference  manuals.  The  other  type  of  error  message  is  generated 
from  within  the  ATB  program,  in  which  case  a  stop  number  will  be  displayed  on  the  computer 
monitor.  You  can  find  the  cause  of  this  stop  by  referencing  the  numbered  stop  list  attached  at  the 
end  of  the  ATB  Input  Manual. 

To  troubleshoot  a  simulation  problem,  the  main  output  .aou  file  often  offers  a  good  indication.  This 
file  first  prints  all  the  input  data  with  labels  for  each  variable.  If  an  error  occurs  in  an  input  card,  the 
output  stops  near  the  error  card.  Therefore,  it  indirectly  indicates  that  there  are  problems  in  the  next 
few  input  cards,  such  as  wrong  format  or  inconsistency  with  data  in  earlier  input  cards. 

Once  the  program  is  ranning,  extensive  dynamic  and  kinematic  data  are  printed  to  the  .aou  file  at 
each  successful  time  step.  These  include  segment  linear  and  angular  positions,  velocities,  and 
accelerations,  joint  forces  and  torques,  and  external  forces  and  torques.  Additionally,  information 
about  when  the  belt  reference  points  are  being  dropped  and  picked  up,  and  convergence  of  the 
integrator  is  also  included  in  the  .aou  file.  If  an  error  occurs,  the  output  stops  and  an  error  message 
is  usually  printed  out  at  the  end  of  the  file. 
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2.  GENERAL  FORMULATION  OF  THE  ATB  MODEL 


The  Articulated  Total  Body  (ATB)  Model  is  primarily  designed  to  evaluate  the  three-dimensional 
dynamic  response  of  a  system  of  bodies  when  subjected  to  a  dynamic  environment  consisting  of 
applied  forces  and  interactive  contact  forces.  Although  the  ATB  Model  was  originally  developed 
to  model  the  dynamic  response  of  crash  dummies  and,  with  later  modifications,  the  response  of  the 
human,  the  ATB  Model  is  quite  general  in  nature  and  can  be  used  to  simulate  a  wide  range  of 
physical  problems  that  can  be  approximated  as  a  system  of  connected  or  free  bodies.  The  ATB 
Model  has  been  used  to  model  such  widely  diverse  physical  phenomena  as  human  body  dynamics, 
the  motion  of  the  balls  in  a  billiards  game,  and  the  transient  response  of  an  MX  missile  suspended 
from  cables  in  a  wind  tunnel.  Version  V  further  expands  the  model’s  capabilities  to  include  water 
force  simulation  with  personal  flotation  devices,  robotic  motion  simulation,  and  deformable  segment 
modeling.  This  flexibility  in  the  ATB  Model  can  cause  the  application  of  the  ATB  program  to 
appear  to  be  overly  complex  to  the  uninitiated  user.  The  purpose  of  this  section  is  to  present  the 
primary  program  features  that  should  be  mastered  to  utilize  Ae  ATB  program.  Throughout  this 
discussion  a  nxmiber  of  input  variables  will  be  mentioned.  A  complete  description  of  these  and  all 
input  variables  is  presented  in  the  ATB  Model  Input  Manual  file  included  with  the  program.  For  a 
more  detailed  discussion  of  how  the  model  uses  these  data,  the  various  theory  manuals  are 
recommended  (1, 3, 4, 9, 10, 11, 18, 13, 14). 

2.1  Chain  Structure  of  the  ATB  Model 

The  system  to  be  simulated  by  the  ATB  Model  can  be  made  up  of  one  or  more  segments  which  may 
be  coimected  or  free.  To  avoid  confusion  between  the  overall  body  or  object  to  be  modeled  and  the 
individual  rigid  or  deformable  bodies  that  make  up  the  overall  body,  throughout  this  report  the  term 
“segment”  will  henceforth  be  used  to  refer  to  the  individual  rigid  or  deformable  bodies  and  the  term 
“body”  will  refer  to  the  overall  body  or  object  to  be  modeled.  The  approach  used  in  the  ATB  Model 
is  to  consider  the  body  as  being  segmented  into  individual  rigid  or  deformable  segments.  Segments 
are  assigned  mass  and  moments-of-inertia  and  are  joined  at  locations  representing  the  physical  joints 
of  the  human  body,  such  as  the  shoulder  joint  or  the  knee  joint. 

The  system  can  be  made  up  of  a  number  of  free  segments,  bodies  of  segments  coupled  together  at 
joints,  or  a  combination  of  both.  A  body  made  up  of  coupled  segments  should  form  an  open  chain 
or  a  tree  structure.  While  this  is  not  an  absolute  requirement,  closed  chains  may  encounter 
computational  problems.  One  must  also  be  careful  not  to  exceed  the  maximum  number  of  segments 
(M/^EG)  specified  by  the  dimension  statements  of  the  program  variables  (See  ATB  Model  Input 
Manual  for  a  list  of  program  parameters  and  their  current  values). 

The  total  number  of  body  segments  (NSEG)  and  total  number  of  joints  (NJNT)  used  to  compose 
all  the  bodies  in  a  simulation  are  input  parameters.  Figure  2  depicts  a  typical  17-segment  (NSEG 
=  17)  model  with  16  joints  (NJNT  =16)  that  is  commonly  used  in  car  crash  and  aircraft  ejection 
simulations  for  modeling  humans  and  dummies.  The  number  of  segments  can  be  readily  varied  in 
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the  input  without  code  modifications  for  up  to  MAXSEG  segments.  Variations  of  the  17- 
segment  body  have  combined  the  forearms  and  hands  for  a  15-segment  body  (NSEG  =15  and  NJNT 
=  14),  added  shoulder  segments,  or  used  different  torso  configurations. 

Whatever  the  specific  body  model,  the  procedure  to  constmct  the  body  remains  the  same.  The  body 
is  assembled  as  a  chain  of  individual  segments.  The  body  can  take  on  a  tree-like  stracture,  with 
several  chains  (here  representing  the  arms  and  the  legs)  branching  out  from  several  connected 
segments. 

The  body  segments  and  joints  are  assigned  identification  numbers,  I  for  the  segments  ranging  from 
1  to  NSEG  and  J  for  the  joints  ranging  from  1  to  NJNT.  The  assignment  of  the  identification 
numbers  is  defined  by  the  order  in  which  the  segments  and  joints  are  listed  in  the  input.  They  are 


used,  along  with  the  one-dimensional  array,  JNT(J)  for  J  =  1  to  NJNT,  to  define  the  connectivity  of 
the  segments.  Segment  1  and  the  first  segment  of  any  body  are  considered  reference  segments. 
Although  the  reference  segment  may  be  any  of  the  body  segments,  it  has  been  found  that,  for  this  17- 
segment  body  model,  the  lower  torso  is  the  best  choice  for  the  reference  segment.  When  the  ATB 
Model  was  first  developed,  the  head  was  chosen  as  the  reference  segment.  It  was  found  that  the 
erratic  accelerations  of  the  head  caused  numerical  problems  with  the  program  integrator  and  that  it 
was  more  beneficial  to  use  a  more  stable,  nonextremity  segment,  hence  the  choice  of  the  lower  torso 
as  the  reference  segment.  A  generalization  of  this  result  is  the  recommendation  that,  regardless  of 
the  body  model,  the  reference  segment  be  chosen  to  be  one  that  undergoes  the  smallest  accelerations 
of  any  of  the  segments  and/or  is  the  heaviest  segment.  Use  of  the  lower  torso  as  the  reference 
segment  also  makes  positioning  the  body  into  a  seat  easier. 

Once  the  lower  torso  is  selected  as  the  reference  segment  and  is  designated  as  segment  number  1, 
the  remaining  segments  should  be  numbered  in  an  order  moving  away  from  the  reference  segment. 
The  left  side  of  Table  1  shows  the  segment  numbers  for  this  17-segment  body.  Each  segment  can 
also  be  given  a  symbol  name  of  up  to  four  alphanumeric  characters  for  output  labeling  purposes,  as 
shown  in  column  3.  The  connectivity  of  the  segments  is  provided  by  the  joints.  The  JNT(J)  array 
provides  the  relationship  between  the  segments  and  joints.  When  JNT(J)  =  I,  the  joint  J  connects 
segment  I  with  segment  J  +  1.  In  other  words,  JNT(J)  stores  the  segment  number  of  the  proximal 
segment  for  the  Jth  joint.  In  the  context  of  the  ATB  Model,  a  proximal  segment  is  the  segment 
nearest  to  the  reference  segment,  whereas  a  distal  segment  is  the  segment  further  from  the  reference 
segment.  For  example  in  Table  1,  joint  J  =  4  (head  pivot)  connects  segment  J  +  1  =  5  (head)  to 
segment  4  (neck).  Therefore  JNT  for  joint  4  is  4.  Another  example  is  joint  J  =  5  (right  hip)  which 
connects  segment  J  +  1  =  6  (right  upper  leg)  to  the  lower  torso,  segment  1.  Therefore,  JNT  for  joint 
Sis  1. 

Successive  segments  and  joints  are  assigned  with  the  provision  that  each  joint  J  connect  segment  J 
+  1  to  a  previously  assigned  segment.  If  the  J  +  1  segment  is  a  reference  segment  for  an  additional 
body,  JNT(J)  is  set  to  0.  This  signifies  that  joint  J  will  be  a  null 
joint  and  that  segment  J  +  1  will  be  the  reference  or  base  segment 
of  another  body.  This  permits  the  specification  of  multiple  bodies 
that  are  disconnected  or  free. 

Besides  using  a  joint,  tWo  segments  can  also  be  connected  by  using 
a  spring-damper  combination,  as  shown  in  Figure  3.  One  situation 
where  you  might  like  to  use  a  spring-damper  combination  is  when 
you  want  to  model  the  thorax  as  two  segments  (spine  and  sternum) 
connected  by  a  spring-damper  combination  in  order  to  evaluate 
chest  deflection. 

Ellipsoids  are  used  to  represent  the  physical  appearance  of  the 
segments.  They  are  the  outer  surfaces  of  the  segments,  and  can 
interact  with  the  environment.  The  B.2  cards  in  the  input  file  are 
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Table  1  Segment  and  Joint  Assignments  and  Connectivity 


■ 

SEGMENT 

NAME 

SYMBOL 

■ 

JOINT 

NAME 

SYMBOL 

JNT(J) 

CONNECTS 

SEGMENTS 

1 

Lower  Torso 

LT 

■; 

1 

Pelvis 

P 

1 

1-2 

2 

Center  Torso 

CT 

2 

Waist 

W 

2 

2-3 

3 

Upper  Torso 

UT 

3 

Neck  Pivot 

NP 

3 

3-4 

D 

Neck 

N 

4 

Head  Pivot 

HP 

4 

4-5 

Head 

H 

5 

Right  Hip 

RH 

1 

1-6 

6 

Right  Upper  Leg 

RUL 

6 

Right  Knee 

RK 

6 

6-7 

7 

Right  Lower  Leg 

RLL 

7 

Right  Ankle 

RA 

7 

7-8 

8 

Right  Foot 

RF 

8 

Left  Hip 

LH 

1 

1-9 

9 

Left  Upper  Leg 

LUL 

9 

Left  Knee 

LK 

9 

9-10 

10 

Left  Lower  Leg 

LLL 

10 

Left  Ankle 

LA 

10 

10-11 

11 

Left  Foot 

LF 

11 

Right  Shoulder 

RS 

3 

3-12 

12 

Right  Upper 

Arm 

RUA 

12 

Right  Elbow 

RE 

12 

12-13 

13 

Right  Lower 

Arm 

RLA 

13 

Right  Wrist 

RW 

13 

13-14 

14 

Right  Hand 

RH 

14 

Left  Shoulder 

LS 

3 

3-15 

15 

Left  Upper  Arm 

LUA 

15 

Left  Elbow 

LE 

15 

15-16 

16 

Left  Lower  Arm 

LLA 

16 

Left  Wrist 

LW 

16 

16-17 

17 

Left  Hand 

LH 

used  to  define  the  segments.  Each  card  contains  the  segment’s  weight,  principal  moments  of  inertia, 
its  ellipsoid’s  semi-axes  and  center  offset  from  the  center  of  gravity  (CG).  The  ellipsoid’s  semi-axes 
define  the  size  of  the  segment.  Following  the  B.2  cards,  a  set  of  B.3  cards  gives  the  location  of  each 
joint.  For  the  Jth  joint,  the  B.3  card  has  the  value  of  JNT(J)  and  the  joint  location  with  respect  to 
each  of  the  two  connected  segments.  The  B  cards  can  be  automatically  generated  using  the  GEBOD 
program  for  human  and  dummy  subjects. 

Although  human  and  dummy  subjects  are  most  often  simulated  in  the  ATB  Model,  there  is  no 
limitation  on  the  type  of  subject.  For  example,  a  MERLIN  robot  arm  modeled  in  six  segments  and 
five  joints  (13)  has  been  simulated.  In  this  case,  the  user  defined  the  chain  stracture  and  built  B.2 
and  B.3  cards  manually  instead  of  using  the  GEBOD  program.  The  program  also  allows  multiple 
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bodies  as  long  as  the  total  segment  number  does  not 
exceed  MAXSEG.  For  example,  two  (driver  and 
passenger)  occupants’  motions  in  an  accident  can  be 
simulated  together.  The  user  needs  only  to  generate  the 
B  cards  for  both  occupants  and  combine  them,  with  a 
null  joint  between.  It  should  be  noted  that  there  are  two 
reference  segments  in  this  case,  the  lower  torso 
segments  of  the  driver  and  the  passenger.  A  body  may 
also  consist  of  only  one  segment.  An  example  is  a 
vehicle  simulation  of  a  pickup  truck’s  rollover  (18). 

The  body  is  the  pickup  truck  represented  by  a  single  , 

segment.  Its  outer  shape  is  described  using  multiple  Figure  4.  One-Segment  Vehicle 

contact  hyperellipsoids  (see  section  2.3.2)  attached  to 
the  segment  (Figure  4)  instead  of  only  one  ellipsoid. 

hi  addition  to  rigid  segments,  the  ATB  Model  accepts  deformable  segments  created  by  finite  element 
analysis  tools.  The  need  for  a  deformable  segment  arises  when  an  accurate  response  is  required  for 
relatively  “flexible”  segments,  such  as  the  human  neck,  or  where  local  segment  vibrations  occur. 
To  use  this  option,  the  user  must  first  develop  a  finite  element  model  of  each  deformable  segment 
and  perform  modal  analysis  to  determine  their  selected  natural  frequencies  and  mode  shapes.  The 
data  required  by  the  ATB  Model  for  each  deformable  segment  are  then  placed  in  separate  data  files. 
Each  deformable  segment’ s  finite  element  data  file  consists  of  node  numbers  and  coordinates,  natural 
frequencies,  and  mode  shapes.  Names  of  these  files  are  input  on  the  B.l.b  card.  The  procedure  for 
developing  deformable  segment  models  is  outlined  in  Appendix  B. 

2.2  Reference  Coordinate  Systems 

The  ATB  Model  utilizes  many  reference  coordinate  systems,  with  respect  to  which  points  in  space 
and  directions  are  calculated  within  the  program.  Considerable  flexibility  in  the  choice  of  coordinate 
systems  and  their  specification  for  both  input  and  output  are  available.  The  primary  coordinate 
systems  used  in  the  model  are  the  inertial,  vehicle,  local  body  segment,  segment  principal  moment 
of  inertia,  joint,  and  contact  ellipsoid  reference  coordinate  systems.  The  specification  of  each 
reference  coordinate  system  requires  an  origin  and  a  direction  cosine  matrix  which  relates  one 
reference  coordinate  system  to  another.  The  direction  cosine  matrix  is  usually  initially  specified  by 
three  rotation  angles,  yaw,  pitch,  and  roll,  as  depicted  in  Figure  5.  These  are  consecutive  body  fixed 
rotations  about  the  Z,  Y,  and  X  axes,  respectively.  All  coordinate  systems  discussed  in  this  section 
are  orthogonal. 
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Figure  5.  Yaw,  Pitch,  and  Roll 


2.2.1  Inertial  Reference  Coordinate  System 

The  ATB  Model  assumes  that  the  coordinates  of  the  origin  of  the  inertial  reference  system  are  zero 
and  all  other  coordinate  systems  are  specified  with  respect  to  this  system.  The  user  may  place  the 
origin  of  the  inertial  reference  coordinate  system  at  any  convenient  point  from  which  his  data  are 
referenced.  The  orientation  is  partially  specified  by  defining  which  way  is  down  by  the  values 
supplied  for  the  components  of  the  gravity  vector.  It  has  been  customary  to  supply  (zero,  zero,  g) 
as  the  components  of  the  gravity  vector  to  specify  that  the  positive  Z,  axis  is  pointing  downward,  as 
shown  in  Figure  6.  Hence,  in  terms  of  a  standing  man,  the  force  of  gravity  would  be  pointing  in  the 
direction  from  his  head  to  his  feet.  Typically,  the  forward  direction  (pointing  from  the  back  of  the 
standing  man  to  his  chest)  is  taken  as  the  positive  Xj  axis  and  (by  the  right  hand  rule)  the  positive 
Yj  axis  is  in  the  lateral  direction  (pointing  from  the  standing  man’s  left  side  to  his  right  side). 


Figure  6.  Inertial  (Ground)  Coordinate  System 
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However,  the  user  may  specify  any  frame  of  reference  that  suits  their  application,  one  with  which 
they  are  more  familiar,  or  in  which  their  input  data  have  been  measured. 

It  is  sometimes  necessary  that  contact  surfaces  (planes  or  ellipsoids)  be  located  with  respect  to  the 
inertial  reference  coordinate  system,  e.g.,  the  ground  for  pedestrian  simulations.  Since  the  program 
assumes  that  contact  surfaces  are  associated  with  segments,  a  special  segment  identification  number 
(NGRND)  is  used  within  the  program  for  this  purpose.  NGRND  is  the  largest  segment  number  used 
by  the  program  and  is  assigned  the  value  NGRND  =  NSEG  +  number  of  airbags  (NBAG)  + 
number  of  vehicles  +  1  and  corresponds  to  the  inertial  coordinate  system.  The  linear  position  and 
velocity  for  this  segment  are  set  to  zero  and  its  direction  cosine  matrix  is  an  identity  matrix 
throughout  the  duration  of  the  simulation.  This  permits  the  use  of  segment  NGRND  for  the 
attachment  of  contact  surfaces. 

2.2.2  Vehicle  Reference  Coordinate  Systems 

Up  to  six  prescribed  motions  can  be  defined.  These  can  be  prescribed  motions  of  segments  defined 
earlier  in  the  B  cards,  or  of  vehicles.  The  primary  vehicle  is  the  last  prescribed  motion  defined  and 
is  different  from  the  other  vehicles  in  that  it  serves  as  the  default  reference  coordinate  system  for 
several  types  of  input  and  output.  Most  of  the  contact  panels  are  usually  defined  with  respect  to  this 
system  and  much  of  the  output  can  be  produced  with  respect  to  this  system.  The  origin  of  each 
vehicle  coordinate  system  is  arbitrary,  and  any  convenient  reference  point  may  be  chosen  to  which 
input  and  output  data  would  be  most  meaningful.  The  frames  of  reference  (the  directions  of  the 
positive  X,  Y,  and  Z  axes)  are  arbitrary  and  should  be  chosen  to  accommodate  input  data. 

A  special  segment  identification  number  is  assigned  for  each  of  the  vehicles  where  NVEHl  =  NSEG 
+1,  NVEH2  =  NSEG  +  2,  etc.  so  that  each  vehicle  may  be  treated  like  other  segments  for  contact 
surface  specifications.  However,  no  matter  how  large  the  computed  contact  forces  and  torques  are 
on  these  vehicle  segments,  the  prescribed  motion  of  the  vehicle  segment  will  not  change.  See  section 
2.5  for  more  information. 

2.2.3  Body  Segment  Reference  Coordinate  Systems 

Each  body  segment  has  a  local  reference  coordinate  system,  sometimes  referred  to  as  the  segment 
geometric  coordinate  system.  Each  body  segment  has  a  mass  and  principal  moments  of  inertia.  The 
local  reference  coordinate  system,  marked  with  subscript  ‘L’  in  Figure  7,  has  its  origin  at  the  segment 
mass  center.  The  principal  moment  of  inertia  axes,  subscripted  with  ‘P’,  are  specified  with  respect 
to  the  local  reference  system.  The  contact  (hyper)  ellipsoid’s  origin  and  orientation,  represented  by 
the  ellipsoid  coordinate  system  in  Figure  8,  are  also  specified  with  respect  to  the  local  reference 
system.  There  is  no  direct  association  within  the  ATB  Model  of  the  segment  inertial  properties  and 
the  contact  (hyper)  ellipsoid.  Unlike  the  vehicle  segments,  a  body  segment’s  kinematics  are 
computed  based  on  the  dynamic  interactions  the  body  segment  experiences  during  a  simulation.  A 
body  segment  can  be  given  an  initial  position,  orientation,  and  linear  and  angular  velocities,  and  its 
motion  is  then  computed  for  the  remainder  of  the  simulation  subject  to  any  imposed  constraints  (e.g. 
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Figure  7.  Segment  Local  and  Principal  Moment  of  Inertia  Coordinate  Systems 

number  and  type  of  joint)  and  applied  forces.  The  motion  of  the  body  segments  cannot  be  specified 
unless  the  body  segment  is  a  reference  segment.  The  orientation  of  the  segment  local  reference 
coordinate  systems  can  be  arbitrarily  defined.  The  standard  convention  has  been  to  choose  the  axes 
so  that  when  the  body  is  in  an  upright  standing  position  with  arms  at  the  side,  the  Z  axis  is 
downward,  the  X  axis  is  to  the  front,  and  the  Y  axis  is  to  the  body’s  right. 

The  segment  center  of  gravity  is  indirectly  determined  through  the  segment’s  joint  coordinates 
because  these  coordinates  are  given  in  the  local  reference  systems  by  Card  B.3.  If  the  segment  is 
deformable,  there  is  no  meaningful  local  reference  system  because  the  center  of  mass  keeps  changing 
when  deformation  occurs.  Its  joint  locations  are  specified  by  the  node  nuiribers. 


12 


A  contact  ellipsoid  is  associated  with  each  segment  and  is  used  for  interactions  with  the 
environment.  The  ellipsoid  number  is  the  same  as  the  segment  number  of  the  associated  segment. 
The  ellipsoid  coordinate  system  origin  is  located  at  the  ellipsoid  center  and  its  axes  are  the  ellipsoid 
semiaxes.  The  ellipsoid  coordinate  system  can  be  translated  from  the  segment  local  coordinate 
system,  but  cannot  be  rotated  in  the  B  cards.  In  order  to  rotate  the  ellipsoid  axes,  the  ellipsoid  must 
be  redefined  in  the  D.5  cards  (additional  (hyper)ellipsoids),  using  the  original  ellipsoid  number.  See 
Section  2.3.2. 

The  dynamic  equations  in  the  ATB  Model  are  solved  in  terms  of  the  principal  axes.  All  three- 
dimensional  bodies  have  an  inertia  tensor.  Six  of  the  nine  inertia  tensor  elements  are  independent, 
therefore  it  is  a  second  order,  symmetric  tensor.  Any  body  has  three  principal  directions  for  which 
there  are  three  moments  of  inertia,  corresponding  to  the  diagonal  elements  of  an  inertia  tensor  when 
all  the  off-diagonal  terms  are  equal  to  zero.  The  segment  principal  coordinate  system  axes 
correspond  to  the  three  principal  directions,  therefore  only  the  three  principal  moments  of  inertia 
must  be  specified. 

The  principal  axes  are  fixed  with  respect  to  the  segment  local  reference  axis.  After  input,  the  ATB 
model  converts  all  data  points  expressed  in  the  local  segment  reference  coordinate  systems  to 
principal  coordinates  and,  prior  to  output,  back  to  the  local  segment  reference  coordinate  systems 
in  a  manner  that  is  transparent  to  the  user.  Therefore,  when  the  input  description  refers  to  local 
segment  reference,  the  local  and  not  the  principal  moment  of  inertia  reference  coordinate  system  is 
implied.  Note  that,  for  some  cases  where  the  principal  axes  are  aligned  with  the  local  reference  axes, 
the  two  are  coincident. 

2.2.4  Joint  Reference  Coordinate  Systems 

In  the  ATB  Model,  the  maximum  number  of  joints  is  MAXJNT.  A  complete  definition  of  a  joint 
consists  of  geometric  location,  joint  coordinate  systems,  type  of  joint,  and  mechanical  properties. 

Based  on  the  mathematical  formulations  in  the  joint  force  and  torque  computation  subroutines,  it  is 
necessary  to  define  two  coordinate  systems  for  a  joint,  one  rigidly  attached  to  each  of  the  two 
segments  that  are  connected  by  the  joint,  as  shown  in  Figure  9.  As  described  in  the  previous 
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sections,  these  two  segments  are  identified  as  segments  JNT  (J)  and  J+1  for  joint  J.  The  origin  of 
each  joint  reference  coordinate  system  (or  the  location  of  the  joint)  is  specified  in  the  segment  local 
reference  coordinate  systems  of  both  segments  JNT  (J)  and  J+1.  The  orientations  of  the  joint  axis 
systems  are  specified  by  rotation  angles  (yaw,  pitch,  and  roll)  from  the  local  reference  systems  of 
both  segments.  Note  that  once  the  two  joint  coordinate  systems  are  defined,  they  are  fixed  in  their 
corresponding  segments  and  do  not  move  relative  to  the  segments.  In  example.ain.  Cards  B.3.a  and 
B.3.b  specify  the  joint  origin  locations  and  joint  coordinate  system  orientations,  respectively. 

In  Card  B.3.a,  the  parameter  IPIN  is  used  to  specify  the  joint  type.  Figure  10  shows  all  the  joint 
types  used  in  the  ATB  Model.  Among  them,  the  null  joint  can  be  used  to  disconnect  two  bodies 
within  the  required  chain  structure.  For  example,  in  an  automobile  crash  simulation  with  two 
occupants,  the  first  occupant’s  joint  data  are  followed  by  a  null  joint  and  then  the  second  occupant’s 
joint  data.  The  ball-and-socket  joint  is  suitable  for  modeling  a  human  shoulder  joint  and  a  pin  joint 
is  suitable  for  an  elbow  joint.  The  Euler  joint  is  a  type  of  joint  which  has  full  three-dimensional 
motion  and  at  the  same  time  allows  the  user  to  impose  various  constraints  on  its  motion.  An  Euler 
joint  is  used  in  the  modeling  of  the  Hybrid  m  dummy’s  shoulder  joint.  The  slip-type  joints  can  be 
used  in  spine  and  neck  modeling,  to  account  for  compression  and  tension. 


Null 


JOINT  TYPES 


Ball  &  Socket  or  Free 


Euler 
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Joint  forces  and  torques  are  computed  by  the  ATB  program  as  a  function  of  the  relative  orientation  of 
the  two  joint  coordinate  systems  at  the  joint,  i.e.,  the  joint  angle  and  angular  velocity.  The  joint 
coordinate  system  associated  with  the  JNT  (J)  segment  is  used  as  the  base  reference  system  for 
determining  the  joint  parameters.  For  the  pin  joint,  the  Y  axis  is  the  axis  of  rotation.  For  the  ball  joint 
or  free  joint,  flexure  (theta)  is  the  angle  between  the  two  Z  axes,  while  azimuth  (phi)  is  the  angle 
between  the  base  X  axis  and  the  projection  of  the  segment  J+l’s  Z  axis  into  the  X-Y  base  plane,  and 
twist  is  rotation  about  the  base  Z  axis.  For  the  Euler  joints,  precession,  nutation  and  spin  are  defined 
as  the  rotations  about  the  Z,  X,  and  Z  axes  respectively  from  the  base  joint  coordinate  system  to  the 
segment  J+l's  joint  coordinate  system.  For  the  slip  joint,  the  linear  motion  is  defined  along  the  base 
Z  axis.  Further  descriptions  of  the  joint  types  and  Aeir  axis  systems  can  be  found  in  References  5  and 
10. 

The  joint  mechanical  properties  define  the  relationship  between  the  joint  resistive  torque  and  the  joint 
angle  and  angular  velocity.  These  properties  include  the  joint  spring,  viscous,  and  coulomb  torque 
characteristics.  Two  options  are  available  for  specifying  the  spring  torques.  Figure  1 1  shows  the  joint 
spring  torque  function  for  the  first  option.  In  this  definition,  a  linear  torque  vs.  joint  angle  is  prescribed 
until  a  specific,  user-defined,  joint  stop  angle,  S5,  is  reached.  For  angles  greater  than  the  joint  stop, 
a  quadratic  or  cubic  restoring  torque  is  added.  Using  this  option,  the  torques  are  symmetric  about  the 
zero  position.  Alternatively,  the  user  has  the  option  to  define  a  joint  resistive  torque  function  that 
depends  on  two  joint  angles,  flexure  and  azimuth,  by  using  the  E.7  card  to  construct  a  two-dimensional 
matrix  data  array.  The  joint’s  B.4  card  then  references  the  function  in  the  E.7  card  for  the  joint  spring 
characteristics  instead  of  using  the  coefficients  in  the  B.4  card.  The  spring  characteristics  of  the  hip 
and  shoulder  joints,  etc.,  in  example.ain  are  defined  using  E.7  cards.  Figure  12  shows  the  viscous  and 
coulomb  torques,  defined  as  functions  of  joint  angular  velocity.  These  functions  are  used  with  both 
spring  torque  options  and  are  defined  in  example.ain  using  the  B.5  cards. 
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Figure  12.  Joint  Torques 


Besides  resistive-type  torque,  the  joint  actuator  option  allows  active  torques  to  be  applied  at  the  joints, 
driving  the  joints  to  prescribed  target  angles.  It  can  be  used  to  simulate  the  joint  actuators  on  a  robot 
or  active  muscles.  The  basic  computational  scheme  uses  feedback  control  logic.  The  F.  10  card  is  used 
to  define  a  relationship  between  the  joint  target  angle  functions  and  their  corresponding  active  torques’ 
control  gain  functions.  For  detailed  information.  Reference  13  is  recommended. 

2.3  Environment  Modeling 

There  are  several  aspects  in  building  the  environment  to  which  the  body  is  exposed.  Contact  planes 
and  ellipsoids  are  used  to  represent  important  geometric  objects  with  which  the  body  may  interact. 
Harness  belts,  airbag  systems,  and  personal  flotation  devices  offer  different  restraint  systems  for  the 
body.  A  predefined  force  or  torque  on  the  body  can  also  be  modeled. 

2.3.1  Modeling  the  Environment  Using  Contact  Planes 

The  contact  planes  are  parallelograms  which  do  not  have  any  inertial  properties,  but  which  provide 
contact  surfaces  used  to  define  the  environment  configuration.  The  interaction  between  the  body  and 
surrounding  environment  can  be  provided  by  contact  forces  between  segments  and  contact  planes. 
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In  the  ATB  Model,  a  set  of  D.2  cards  is  used  to  define  the  contact  planes.  A  plane  is  represented  by 
the  coordinates  of  its  three  comer  points,  Pj,  P2,  and  P3  shown  in  Figure  13.  These  coordinates  are 
given  in  the  reference  system  of  the  segment,  vehicle  segment,  or  the  ground  (inertial  segment)  to 
which  the  plane  is  attached.  The  segment  to  which  the  plane  is  attached  is  specified  in  the  F.l  cards. 
See  Section  2.4  for  further  explanation.  The  order  in  which  these  points  are  listed  defines  the  positive 
side  of  the  plane.  If  the  order  is  Pj,  Pj,  and  P3,  then  the  force  generated  by  the  plane  contact  will  be 
in  the  direction  of  the  vector  produced  by  the  cross-product  PjPj  x  P1P3,  as  shoAvn  in  Figure  14.  Users 
can  refer  to  the  D.2  cards  in  example.ain  for  an  example  of  modeling  a  seat  using  contact  planes. 


A 


Figure  14.  Positive  Side  of  Plane 
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The  maximum  number  of  contact  planes  is  MAXPLN.  Though  a  user  can  use  as  many  planes  as 
possible  to  model  the  environment,  the  number  of  plane  contacts  can  affect  the  computation  time. 

2.3.2  Additional  Contact  (Hyper)Ellipsoids 


In  addition  to  the  ellipsoids  defined  with  the  body  segments,  the  ATB  Model  has  an  option  to  attach 
contact  (hyper)ellipsoids  to  body  segments,  vehicle  segments  or  the  ground  (inertial)  segment.  The 
(hyper)ellipsoids  have  no  mass  or  moments  of  inertia  and  hence  no  dynamic  response.  They  are  rigidly 
attached  to  a  segment  at  a  point  and  with  an  orientation  specified  with  respect  to  the  segment  local 
reference  coordinate  system,  as  shown  in  Figure  15.  As  with  the  planes,  the  segment  to  which  the 
(hyper)ellipsoid  is  attached  is  chosen  on  the  F.l  or  F.3  cards.  See  Section  2.4.  The  contact 
(hyper)ellipsoid  coordinate  system  is  formed  by  the  three  orthogonal  semi-axes  of  the  (hyper)ellipsoid, 
with  the  coordinate  system’s  origin  at  the  geometric  center  of  the  (hyper)ellipsoid.  A  normal  contact 
ellipsoid  has  its  power  values,  m,  n,  and  p,  equal  to  2.  An  ellipsoid  with  its  powers  greater  than  2  is 
called  a  hyperellipsoid.  The  function  describing  a  (hyper)ellipsoid  is  (x/a)®  +  (y/b)°  +  {zlcf  =  1.  As 
the  powers  increase,  the  shape  of  the  hyperellipsoid  becomes  more  square.  An  example  of  using 
contact  hyperellipsoids  is  the  modeling  of  a  pickup  tmck’s  squared  exterior  shape,  as  shown  in  Figure 
4.  Like  the  planes,  the  (hyper)ellipsoids  are  for  contact  purposes  only. 

The  location  of  the  contact  (hyper)ellipsoid  is  specified  with  an  offset  vector  which  starts  at  the  origin 
of  the  segment’s  local  reference  coordinate  system  and  ends  at  the  point  where  the  center  of  the  contact 
(hyper)ellipsoid  is  to  be  attached.  The  orientation  of  the  contact  (hyper)ellipsoid  is  specified  by 
rotation  angles  with  respect  to  the  local  segment  reference  system.  If  no  rotation  angles  are  specified 
for  the  contact  (hyper)ellipsoid,  the  X,  Y,  and  Z  semi-axes  of  the  contact  (hyper)ellipsoid  will  coincide 
with  the  X,  Y,  and  Z  axes  of  the  local  reference  system  of  the  segment. 


Figure  15.  Additional  (Hyper)Ellipsoids 


Contact  (hyper)ellipsoids  can  be  used  for  (hyper)ellipsoid/(hyper)ellipsoid  contacts  and  (hyper) 
ellipsoid/plane  contacts,  but  only  normal  ellipsoids  can  be  used  for  belt/ellipsoid  contacts,  harness 
belt/ellipsoid  contacts  or  airbag/ellipsoid  contacts  (where  the  airbag  is  a  special  type  of  contact 
ellipsoid).  More  than  one  contact  (hyper)ellipsoid  can  be  attached  to  one  (body,  vehicle,  or  ground) 
segment. 

2.3.3  Belt  Restraint  Systems 

The  ATB  Model  provides  two  options  for  modeling  of  belt  restraint  systems:  simple  belt  and  harness 
belt  systems.  In  a  simple  belt  system,  each  restraint  belt  is  assumed  to  lie  in  a  plane  defined  by  two 
anchor  points  attached  to  a  segment  (usually  the  vehicle)  and  by  a  fixed  point  on  a  contact  ellipsoid 
rigidly  attached  to  some  other  segment  (see  Figure  16).  Therefore,  the  belt  is  restricted  to  pass  around 
a  single  segment.  Although  several  simple  belts  may  be  used  in  an  application,  they  cannot  interact 
with  each  other.  The  dynamic  properties  of  the  simple  belt  are  defined  by  initial  slack,  a  force-strain 
function,  and  friction  of  the  contact  between  the  belt  and  the  segment’s  ellipsoid.  However,  the 
friction  is  limited  to  either  zero  or  infinite.  A  strain-rate-dependent  function  is  not  allowed  for  the 
simple  belt.  The  main  limitation  of  the  simple  belt  model  is  that  the  point  at  which  the  belt  contacts 
a  segment  is  fixed  to  the  segment  and  moves  with  it.  The  simple  belt  system  is  modeled  by  a  set  of 
D.3  and  F.2  cards  in  input  files. 


Contact 

Point 


Points 

Figure  16.  Simple  Belt 

The  harness  belt  model  overcomes  some  of  the  simple  belt’s  limitations  by  allowing  interactive  belts 
that  can  slip  over  multiple  segments.  A  harness  consists  of  one  or  several  belts.  Each  belt  is  formed 
by  a  set  of  straight  line  segments  connecting  prescribed  reference  points,  shown  in  Figure  17. 
Endpoints  of  the  belt  may  be  anchor  points  or  junction/tie  points  where  several  belts  may  join 
together.  In  example.ain,  the  set  of  F.8  cards  models  a  harness  system  of  four  belts:  a  conventional 
double  shoulder  strap  and  a  negative  G  strap  tied  together  at  the  middle  point  of  a  lap  belt. 
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Belt  Line 


Figure  17.  Harness  Belt 

The  harness  belt  reference  points  are  points  of  contact  between  the  belts  and  contact  ellipsoids.  Their 
X,  Y,  and  Z  coordinates  in  the  contact  ellipsoid  reference  systems  are  given  in  Card  F.S.dl  to 
determine  their  location  on  the  contact  ellipsoids’  surfaces.  These  coordinates  are  the  only  ATB  input 
data  specified  in  terms  of  the  contact  ellipsoid  coordinate  systems.  The  supplied  values  are  adjusted 
by  the  program  to  lie  on  the  ellipsoid  surface.  Additional  contact  ellipsoids  can  be  attached  to  a 
segment  to  better  model  the  surface  of  the  body.  For  example,  in  modeling  the  shoulder  belt  of  a  three- 
point  harness  shown  in  Figure  18,  an  additional  contact  ellipsoid  with  reference  points  is  attached  to 
the  upper  left  portion  of  the  upper  torso  to  represent  the  belt  layout  on  a  human  shoulder.  It  should 
be  pointed  out  that  hyperellipsoids  cannot  contact  harness  belts. 
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Depending  upon  the  layout  of  the  belt  during  each  time  step,  some  reference  points  may  be  “dropped” 
from  the  calculation  of  the  belt  trajectory  and  forces.  A  surface  normal  is  used  to  determine  whether 
a  reference  point  will  be  included  in  the  calculations.  As  shown  in  Figure  17,  the  surface  normal  is 
an  outward  normal  vector  to  the  surface  of  the  ellipsoid  at  the  reference  point.  If  the  net  belt  force  on 
the  ellipsoid  at  this  reference  point  has  a  positive  component  along  this  normal,  the  point  will  be 
ignored  in  computing  the  belt  forces.  These  dropped  points  may  be  picked  up  at  a  later  time.  In 
addition,  if  no  ellipsoid  is  specified  for  a  reference  point,  this  point  will  always  be  used  in  the 
calculations.  The  user  can  find  those  reference  points  being  dropped  and  picked  up  at  any  particular 
time  step  by  checking  the  primary  output  file  *.aou,  which  is  example.aou  in  the  example  case.  A 
detailed  description  of  the  belt  algorithms  is  given  in  Reference  3. 

For  a  harness  belt,  you  are  required  to  define  a  strain  or  strain-rate-dependent  force  function  for 
computing  belt  forces,  ff  the  belt’s  reference  points  are  allowed  to  penetrate  into  the  surface  of  the 
ellipsoids,  force-deflection  functions  are  also  needed  to  describe  these  penetrations.  These  functions 
are  assigned  to  each  reference  point  in  Cards  F.8.d.  If  no  force-deflection  functions  are  provided,  the 
ATB  Model  assumes  the  surface  is  rigid  and  no  perturbations  of  the  reference  points  normal  to  the 
surface  are  allowed.  Furthermore,  users  can  specify  initial  slack  where  a  negative  value  indicates  a 
pre-tightened  belt. 

A  reference  point  can  move  tangential  to  the  ellipsoid  surface,  both  along  and  normal  to  the  belt  line. 
Friction  coefficients  between  the  belt  and  each  ellipsoid  can  be  defined  in  order  to  control  the  belt 
movement  on  the  surface.  Once  the  belt  contact  ceases  at  a  reference  point,  the  reference  point  will 
remain  at  its  last  belt  contact  position  on  the  surface,  until  it  has  a  negative  normal  force  and  is  picked 
up  again.  If  the  belt  separates  completely  firom  an  ellipsoid  and  contacts  it  again  later,  the  belt  may  not 
be  able  to  pick  up  any  reference  points  if  the  ellipsoid  has  rotated  significantly.  In  that  case,  the  belt 
may  cut  through  the  ellipsoid  without  producing  any  resistance  forces.  Therefore,  in  a  simulation 
involving  complex  body  motion,  some  additional  reference  points,  such  as  the  dropped  point  in  Figure 
17,  may  need  to  be  defined  in  the  input  file.  These  will  be  used  solely  for  later  contacts.  Care  must 
be  taken  to  make  sure  these  points  will  not  be  picked  up  at  the  beginning  of  the  simulation  in  order  to 
avoid  an  unrealistic  belt  configuration. 

2.3.4  Simple  Airbag  Restraint  System 

The  simple  airbag  model  is  a  non-stretchable  bag  of  ellipsoidal  shape  which  interacts  with  contact 
ellipsoids  attached  to  selected  segments  of  the  bodies  or  the  vehicle.  Figure  19  gives  a  complete 
picture  of  the  airbag  model.  Those  contact  ellipsoids  attached  to  the  vehicle  for  holding  up  or  confining 
the  airbag  are  called  reaction  panels.  At  least  one  such  contact  ellipsoid,  called  the  primary  reaction 
panel,  is  required  in  modeling  the  airbag.  A  point  on  this  panel  is  specified  as  the  deployment  point 
from  which  the  bag  deploys.  At  the  beginning  of  the  simulation,  the  bag  is  assumed  to  have  zero 
volume  and  be  located  at  the  deployment  point.  After  a  specified  time,  the  bag  is  inflated  by  using  the 
gas  dynamic  relations  for  the  choked  flow  of  gas  through  a  nozzle.  The  gas  source  is  a  high  pressure 
tank  of  constant  volume.  The  total  amount  of  gas  coming  through  the  nozzle  is  the  volume  of  gas  in 
the  fully  inflated  bag,  at  atmospheric  pressure. 
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During  inflation,  the  size  of  the  bag  is  determined  by  scaling  the  semi-axes  of  the  ellipsoid  by  the  cube 
root  of  the  volume.  The  center  of  the  bag  lies  on  a  vector  which  has  one  end  at  the  deployment  point 
and  is  parallel  to  the  X  axis  of  the  primary  reaction  panel  but  in  the  negative  X  direction.  The  center’s 
distance  from  the  deployment  point  is  equal  to  the  X  semi-axis  of  the  sealed  bag  at  that  instant  of  time. 
The  set  of  ellipsoids  drawn  by  dashed  lines  in  Figure  19  illustrates  this  inflation  process. 

When  the  bag  is  fully  inflated,  it  is  assigned  mass  properties  and  moves  dynamically  like  any  other 
mass  system.  Until  it  is  fully  inflated,  the  orientation  of  the  bag  with  respect  to  the  vehicle  is  held 
constant  and  equal  to  its  initial  orientation.  The  dynamic  motion  of  the  bag  is  updated  by  the  program 
integrator.  An  artificial  spring  force  is  applied  at  the  end  of  the  positive  X  axis  of  the  bag  and  is 
exterior  to  the  primary  reaction  panel.  This  holds  the  bag  to  the  panel. 

During  the  inflation,  the  bag  is  assumed  to  be  at  atmospheric  pressure  and  hence  no  contact  forces  are 
produced  until  the  sum  of  the  instantaneous  volume  of  the  bag  and  the  volume  of  intersection  due  to 
contacts  with  other  segments  reaches  the  geometric  volume  of  the  fully  inflated  bag.  Once  the  bag  is 
fully  inflated,  any  additional  gas  from  the  gas  tank  or  an  increase  in  the  volume  of  intersection  will 
cause  the  pressure  in  the  bag  to  increase  and  thus  produce  contact  forces  on  any  contact  ellipsoids 
intersecting  the  bag.  Each  intersection  of  a  contact  ellipsoid  and  the  bag  is  treated  separately  by  the 
ATB  program,  which  computes  the  decrease  in  the  volume  of  the  bag,  the  effective  area  of  the  contact 
and  the  force  and  torque  per  unit  pressure.  After  all  the  contacts  have  been  considered,  the  total 
decrease  in  volume  is  used  to  compute  the  pressure  of  the  gas  in  the  bag  and  then  the  forces  and 
torques  are  applied  to  the  various  ellipsoids  at  their  maximum  penetration  point  into  the  bag. 

A  set  of  D.4  cards  are  used  in  the  ATB  Model  to  define  airbag  parameters.  A  detailed  airbag 
formulation  can  be  found  in  Reference  5. 
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2.3.5  Applied  Force  and  Torque 


The  ATB  model  has  the  capability  to  apply  time-dependent  forces  and  torques  to  body  segments,  as 
shown  in  Figure  20.  A  force/torque  coordinate  system  is  defined  such  that  a  positive  force  is  applied 
in  the  positive  X  direction  of  the  force/torque  coordinate  system  and  a  positive  torque  is  applied  about 
the  positive  X  axis  of  the  force/torque  coordinate  system  using  the  right-hand  rale.  The  origin  and 
orientation  (rotation)  of  the  force/torque  coordinate  systems  are  specified  with  respect  to  the  local 
reference  coordinate  system  of  the  segment  to  which  the  force/torque  is  to  be  applied.  Cards  D.9 
specify  these  parameters  for  the  ATB  simulation. 


Figure  20.  Applied  Forces 


2.3.6  Wind  Force  Modeling 

The  wind  force  option  was  developed  to  tackle  the  case  where  pilot  ejection  is  simulated.  It  applies 
pressure  type  forces,  such  as  aerodynamic  forces,  to  any  segments  which  penetrate  a  boundary  plane, 
called  the  wind  plane,  as  shown  in  Figure  21 .  Once  a  segment’s  ellipsoid  penetrates  the  wind  plane, 
an  estimate  of  the  projected  area  normal  to  the  wind  pressure  is  made,  and  the  force  and  torque  are 
computed  and  applied  to  the  segment.  For  partial  penetration,  the  force  is  applied  at  the  center  of  the 
intersection  ellipse  between  the  ellipsoid  and  wind  plane.  At  full  penetration,  the  force  is  applied  at 
the  center  of  the  ellipsoid. 

The  wind  plane  is  defined  using  the  D.2  card,  the  same  as  other  regular  planes.  The  user  must  state 
explicitly  in  the  F.7.a  card  which  segments  are  desired  for  wind  force  calculation.  Then,  the  F.7.b  card 
associates  these  segments  with  the  wind  plane  and  wind  pressure  functions,  as  well  as  drag  coefficients 
defined  in  Card  E.6.  There  are  two  types  of  wind  pressure  functions.  The  first  is  a  time-dependent 
wind  pressure  function  which  gives  the  X,  Y  and  Z  components  of  the  wind  pressure  vector  in  the 
reference  coordinate  system  with  respect  to  time.  The  second  type  computes  the  wind  pressure  vector 
as  a  function  of  the  relative  velocity  of  a  segment.  A  time-dependent  drag  coefficient  function  can  also 


23 


^  ^  Wind  Plane 

Figure  21.  Wind  Forces 

be  defined  in  Card  E.6.  The  wind  force  acting  on  a  segment  will  equal  the  wind  pressure  function 
value  multiplied  by  the  drag  coefficient.  If  no  drag  coefficient  is  defined,  1.0  will  be  used. 
Additionally,  a  method  of  calculating  the  wetted  area  to  account  for  segments  blocking  the  wind  is  also 
available.  A  detailed  formulation  of  the  wind  force  modeling  method  can  be  found  in  References  3 
and  10. 

2.3.7  Water  Force  Simulation  Environment  Modeling 

The  water  force  simulation  is  used  to  predict  gross  human  or  dummy  body  response  due  to  water 
forces  with  the  body  fully  or  partially  submerged  and  with  or  without  attached  PFDs  (Personal 
Floatation  Devices).  Modeling  of  water  forces  requires  representation  of  the  water  surface  and  PFDs. 
Figure  22  shows  a  human  subject  wearing  PFDs  floating  above  the  water  surface.  The  subject  is 
partially  submerged. 


Figure  22.  Water  Force  Simulation  of  a  Human  Subject  with  PFDs 
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A  representation  of  the  water  surface  consists  of  the  mean  water  surface  and  waves  forming  the  free¬ 
water  surface.  The  mean  water  surface  is  defined  as  the  X-Y  plane  of  a  Cartesian  coordinate  system 
called  the  water  frame.  The  water  frame’s  origin  and  orientation  are  given  with  respect  to  the  inertial 
coordinate  system,  and  its  Z  axis  points  downwards  into  the  water.  Figure  23  depicts  this  relationship. 
Therefore,  the  mean  water  surface  is  fixed  in  space  and  time  and  used  as  the  location  of  the  water 
surface.  To  have  both  spatial  and  time  variations  for  the  actual  free-water  surface,  waves  are  defined 
and  superimposed  on  the  mean  water  surface.  There  are  two  representations  of  the  free-water  surface 
available  in  the  ATB  Model.  One  uses  a  set  of  regular  waves  which  are  two-dimensional,  sinusoidal 
waves.  The  user  inputs  wave  length,  amplitude  and  phase  angle,  etc.,  to  define  these  waves.  The  ATB 
Model  allows  the  user  to  utilize  up  to  ten  regular  waves  to  describe  the  free-water  surface.  The 
program  will  superimpose  the  components  due  to  each  wave  in  computing  the  free-water  surface. 
Another  option  aJlows  the  user  to  represent  the  free-water  surface  by  a  single  regular  wave  based  on 
the  Pierson-Markovitz  spectrum  for  fully-developed  ocean  waves.  The  user  only  needs  to  supply  a 
wind  velocity  at  a  standard  height  of  63.98  ft  above  the  free  surface.  The  ATB  Model  will  compute 
the  rest  of  the  wave  parameters. 


For  the  portion  of  the  body  submerged  in  water,  the  ATB  Model  computes  the  water  forces  and  torques 
acting  on  it.  The  water  forces  include  the  effects  of  hydrostatic  pressure,  wave  excitation,  added  mass, 
drag  and  lift.  The  hydrostatic  effects  arise  as  a  result  of  hydrostatic  fluid  pressure  acting  on  a  body. 
The  wave  excitation  effects  are  due  to  the  dynamic  pressure  exerted  by  the  waves.  The  added  mass 
effects  result  from  a  volume  of  surrounding  fluid  accelerated  with  a  body.  The  parameters  describing 
these  effects  are  supplied  in  the  ATB  input  file  by  a  set  of  corresponding  coefficients. 

The  ATB  Model  allows  the  user  to  approximate  each  PFD  by  up  to  five  rigid  ellipsoids.  These 
ellipsoids  are  called  PFD  ellipsoids.  Each  PFD  ellipsoid  is  modeled  as  if  it  is  rigidly  connected  to  a 
body  segment  and  has  similar  characteristics  as  a  regular  contact  ellipsoid.  The  semi-axes,  ellipsoid 
center  offset,  and  orientation  of  the  ellipsoid  coordinate  system  are  defined  with  respect  to  the  local 
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reference  system  of  the  segment  to  which  it  is  attached.  The  ATB  Model  allows  modeling  of  up  to  five 
PFDs. 


In  an  ATB  input  file,  a  set  of  F.9  cards  are  used  to  define  the  water  surface,  water  force  coefficients, 
PFDs,  and  other  water  force  simulation  parameters.  For  a  detailed  theoretical  formulation  of  the  water 
force  model.  Reference  12  is  recommended. 

2.4  Contact  Definitions 

The  interactions  between  the  body  and  the  environment  are  expressed  in  terms  of  the  contacts  which 
happen  between  the  body’s  segments  and  the  elements  representing  the  environment.  There  are  seven 
major  types  of  contacts  used  in  the  ATB  Model; 

1.  Plane/ellipsoid 

2.  Ellipsoid/ellipsoid 

3.  Segment/belt 

4.  Segment/airbag 

5.  Segment/water 

6.  Segment/hamess  belt 

7.  Segment/spring-damper 

The  last  five  types  of  contacts  have  been  described  in  previous  sections.  This  section  focuses  on  the 
contact  models  used  most  frequently  (plane/ellipsoid  and  ellipsoid/ellipsoid  contacts)  and  the  functions 
used  by  all  of  the  contact  models. 

2.4.1  Plane/Ellipsoid  and  Ellipsoid/Ellipsoid  Contact 

The  outer  shape  of  each  segment  is  defined  by  the  contact  ellipsoid  attached  to  the  segment.  Most 
segments  have  only  one  contact  ellipsoid  while  some  segments  may  have  several  contact  ellipsoids, 
hyperellipsoids,  or  planes.  Since  the  vehicle  and  ground  are  also  segments,  they  may  have  contact 
(hyper)ellipsoids  attached  to  them,  as  well  as  planes.  Therefore,  a  contact  with  the  segment  is  in  effect 
a  contact  with  the  segment’s  contact  (hyper)ellipsoids.  Referring  to  example.ain,  the  sets  of  Cards  F.l 
and  Cards  F.3  define  the  plane/ellipsoid  and  ellipsoid/ellipsoid  contacts,  respectively,  in  this  way.  The 
user  can  define  contacts  for  all  of  the  combinations  between  planes  and  ellipsoids,  and  ellipsoids  and 
ellipsoids.  However,  it  will  save  computational  time  to  define  only  the  contacts  which  are  likely  to 
occur.  Several  test  mns  can  be  used  to  refine  and  reduce  the  contact  definitions. 
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For  a  plane/ellipsoid  contact,  the  contact  forces  consist  of  a  normal  force  and  a  friction  force  computed 
by  the  ATB  Model’s  force-deflection  routines.  When  an  ellipsoid  contacts  with  a  plane,  it  penetrates 
the  plane.  The  ATB  Model  finds  the  maximum  penetration  of  the  given  ellipsoid  into  the  given  plane 
at  each  time  step,  as  shown  in  Figure  24.  This  penetration  is  the  deflection  meant  by  the  ATB  Model. 
The  contact  forces  corresponding  to  this  deflection  are  calculated  from  the  contact  property  functions 
used  in  Cards  F.3  for  defining  this  contact.  The  computed  forces  are  applied  at  a  point  along  the  line 
joining  the  point  of  maximum  penetration  and  the  center  of  the  intersection  area.  The  ATB  Model  has 
the  capabilities  to  conduct  an  edge  effect  test,  handle  complete  penetration  by  the  ellipsoid,  and  deal 
with  infinite  planes.  The  edge  effect  test  handles  the  situation  where  only  part  of  the  intersecting 
ellipse  is  within  the  plane  boundaries.  The  infinite  plane  option  allows  the  user  to  assume  an  infinite 
plane;  therefore,  no  boundary  test  is  made.  It  should  be  noted  that  the  plane  has  a  positive  side  and 
a  negative  side,  as  described  earlier.  The  contact  force  vector  has  the  same  direction  as  the  plane 
normal,  directed  out  of  the  positive  side.  If  an  ellipsoid  comes  into  contact  from  the  negative  side  of 
the  plane,  it  results  in  a  sudden  large  contact  force  pushing  the  ellipsoid  through  the  positive  side  of 
the  plane. 


Similar  to  plane/ellipsoid  contacts,  the  contact  forces  generated  from  an  ellipsoid/ellipsoid  contact  are 
functions  of  the  penetration  value  of  one  ellipsoid  into  another,  as  shown  in  Figure  25.  The  penetration 
value  is  decided  by  contracting  both  ellipsoids  until  a  single  point  of  contact  is  achieved.  Figure  25 
shows  an  exterior  contact  in  which  one  ellipsoid  approaches  another  from  the  exterior.  The  ATB 
Model  also  allows  an  interior  contact  in  which  ellipsoid  A  contacts  ellipsoid  B  at  B’s  interior  surface. 
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though  this  is  a  rarely  used  option.  For  interior  contacts,  ellipsoid  A  must  be  completely  inside 
ellipsoid  B  before  contact.  _ 


Figure  25.  Ellipsoid-Ellipsoid  Contact 


It  should  be  pointed  out  that  the  contact  properties  are  mutual  characteristics  associated  with  each 
specific  paired  contact.  For  example,  if  test  data  are  available,  the  contact  properties  between  the  head 
ellipsoid  and  dashboard  plane  should  be  different  from  those  between  the  upper  torso  ellipsoid  and 
dashboard  plane. 

2.4.2  Functions  of  Contact  Properties 

For  many  of  the  contact  definitions,  a  set  of  function  numbers  are  used  to  define  the  contact  properties 
for  that  particular  contact.  Each  individual  contact  function  is  defined  in  Cards  E.l  through  E.4. 
Contact  properties  are  described  by  a  combination  of  individual  functions.  In  most  cases  there  are  two 
ways  to  specify  this  combination. 

The  first  method  uses  five  functions  in  combination  to  describe  the  contact  properties.  They  are:  the 
base  force-deflection  function,  the  inertial  spike  function,  the  energy  absorption  function,  the 
permanent  deflection  function,  and  the  firiction  coefficient  function.  The  function  numbers  of  these 
five  functions  are  used  in  the  contact  definition.  Figure  26  depicts  this  concept.  These  are  all 
functions  of  deflection  or  constant  values.  The  base  force-deflection  function  is  used  to  determine  the 
normal  contact  force.  The  inertial  spike  function  is  used  to  model  the  inertial  loading  that  might  take 
place  when  a  plane/ellipsoid  contact  is  initiated.  The  definition  of  each  inertial  spike  function  includes 
an  abscissa  value  D3  in  Card  E.2  for  the  function  such  that  if  unloading  occurs  after  deflection  exceeds 
D3,  the  inertial  spike  is  to  be  ignored.  An  example  using  inertial  spike  might  be  the  contact  between 
the  head  and  the  car’s  windshield,  in  which  case  an  extra  inertial  loading  is  needed  before  the 
windshield  is  broken.  The  energy  absorption  function,  also  called  the  R  factor,  ranges  from  0  to  1  and 
is  used  to  specify  the  amount  of  energy  recovered  at  the  end  of  unloading.  The  permanent  deflection 
function,  also  called  the  G  factor,  ranges  from  0  to  1  and  is  used  to  model  permanent  deformation 
XpERM  due  to  the  contact  force.  For  the  subsequent  contacts,  loading  will  not  start  until  Xperm-  Both 
the  R  and  G  factors  are  used  to  approximate  the  effects  of  hysteresis,  defining  the  unloading  and 
reloading  curve  calculations.  The  unloading  curve  is  a  quadratic  polynomial  from  the  base  curve  to 
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XpERM-  The  reloading  curve  is  a  cubic  pol5aioniial  from  the  point  of  reloading  to  the  base  curve  at 
Xmax  The  friction  coefficient  function  is  used  to  compute  the  contact  friction  force  which  is 
proportional  to  the  normal  force  and  in  the  opposite  direction  of  the  tangential  velocity.  In  summary, 
this  method  establishes  a  contact  behavior  as  first  loading  along  the  base  force-deflection  curve  plus 
the  inertial  spike  (if  it  exists),  then  proceeding  down  an  unloading  curve  between  XpgRM  and  X^ax 
the  deflection  reaches  X^ax-  In  example.ain,  the  ellipsoid/ellipsoid  contact  definitions  in  Cards  F.3 
use  this  method. 


The  second  method  describes  the  contact  properties  using  a  rate-dependent  function,  F(x,  x'),  where 
X  and  x'  are  the  deflection  and  deflection  rate,  respectively.  A  combination  of  four  individual  functions 
are  used  to  construct  F  as  shown  in  Figure  27.  Fj,  F2,  F3,  and  F4  do  not  have  to  have  any  special 
physical  meanings  and  may  be  used  as  pure  mathematical  expressions  for  the  purpose  of  constructing 
F.  The  function  numbers  of  these  four  individual  functions  plus  the  friction  coefficient  function  are 
used  in  the  contact  definition.  This  option  is  invoked  by  setting  the  function  numbers  of  F2,  F3,  and 
F4  in  Cards  F.l  or  F.3  all  to  be  negative.  In  example.ain,  the  block  of  F.l  cards  defining  plane/ellipsoid 
contacts  uses  this  method. 
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F=F,(x')  +  F,(x')F3(x')  +  F,(x') 


Figure  27.  Rate-Dependent  Functions 


Each  individual  function  used  in  the  above  methods  can  be  defined  using  tabular  data  and/or  a 
polynomial  expression  of  up  to  the  fifth-degree  in  Cards  E.l  to  E.4.  Constant  value  functions  are 
allowed.  For  example,  in  example.ain,  function  14  is  a  constant  function  specifying  the  friction 
coefficient;  function  13  is  a  tabular  function  defining  a  base  force-deflection  function  for  contact 
between  ellipsoids  and  stiff  surfaces;  function  31  is  a  linear  (first-order  polynomial)  function  defining 
the  force-strain  relation  for  harness  belts.  Furthermore,  each  individual  function  may  be  subdivided 
into  two  adjacent  functions  fj  and  f2,  where  the  upper  abscissa  value  of  fj  will  be  the  lower  abscissa 
value  of  f2,  as  indicated  in  Figure  28.  The  input  formats  of  fj  and  f2  depend  on  the  signs  of  DO,  Dl, 
and  D2  on  the  E.2  cards. 
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2.5  Prescribed  Motion  and  Initial  Conditions 
2.5.1  Vehicle  Definition  and  Prescribed  Motion 

A  vehicle  in  the  ATB  Model  is  a  massless  segment  with  prescribed  motion.  Unlike  the  body  segments 
described  by  B  cards,  vehicle  segments  are  automatically  created  by  C  cards  where  its  prescribed 
motion  is  defined.  Vehicle  segments  are  primarily  used  to  offer  a  reference  coordinate  system  (see 
section  2.2.2)  and  establish  an  acceleration  time  history  for  the  environment  in  which  the  bodies  reside. 
Therefore,  it  usually  has  contact  planes  and/or  additional  ellipsoids  defined  with  respect  to  it  (see 
section  2.3).  For  example,  in  example.ain,  the  sled  is  a  vehicle  segment  with  a  deceleration  sequence 
defined  in  the  C.3  cards.  The  seat  panels  are  defined  in  the  vehicle  coordinate  system,  i.e.,  with  respect 
to  the  sled. 

The  ATB  Model  allows  up  to  six  prescribed  motions  to  be  defined.  Each  prescribed  motion  has  its 
own  set  of  C  cards,  with  the  primary  vehicle  being  defined  by  the  last  set  of  C  cards  in  the  input  file. 
The  remaining  prescribed  motions  not  associated  with  a  segment  defined  in  the  B  cards  are  called 
secondary  vehicles.  The  numbering  convention  for  segments  in  the  ATB  Model  is,  from  low  to  high, 
the  body  segments  followed  by  the  secondary  vehicles,  the  primary  vehicle,  the  airbags  if  any,  and  then 
the  inertial  system  segment,  which  is  the  ground  segment  by  default.  In  example.ain,  the  sled  is  the 
only  vehicle;  therefore  it  is  also  a  primary  vehicle.  The  secondary  vehicles  are  very  useful  in  modeling 
objects  which  move  within  the  primary  vehicle.  For  example,  in  a  pickup  truck  simulation  with  roof 
crush  as  shown  in  Figure  29,  the  primary  vehicle  is  the  pickup  truck  and  the  roof  is  modeled  as  a 
secondary  vehicle.  Most  planes  are  defined  with  respect  to  the  primary  vehicle,  which  has  its  segment 
local  reference  system  at  the  point  described  by  the  pickup  truck’s  prescribed  motion.  However,  the 
roof  and  side  rails  are  attached  and  defined  in  the  local  reference  system  of  the  secondary  vehicle, 
which  has  the  time  history  of  roof  crush  movement  as  the  prescribed  motion.  The  roof  crush 
movement  in  turn  may  be  defined  with  respect  to  the  primary  vehicle  or  the  inertial  reference  system. 


Figure  29.  Pickup  Truck  Rollover  with  Roof  Crush 
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The  ATB  Model  requires  at  least  one  vehicle,  the  primary  vehicle,  to  be  defined  in  every  simulation. 
If  there  is  no  object  suitable  for  vehicle  definition,  a  motionless  dummy  vehicle  can  be  defined  as  the 
primary  vehicle  for  the  simulation.  The  water  force  simulation  in  Section  2.3.7  is  such  a  case. 

A  body  segment  defined  in  the  B  cards  may  also  be  given  prescribed  motion  if  the  segment  is  a 
reference  segment.  Similar  to  a  vehicle,  the  body  segment’s  prescribed  motion  is  defined  by  a  set  of 
C  cards. 


Each  prescribed  motion  time  history  in  the  C  cards  is  specified  relative  to  a  reference  system.  The 
reference  system  can  be  that  of  another  prescribed  motion  segment,  or  the  default  ground  (inertial) 
segment.  Four  options  are  available  for  specifying  the  motion  data.  They  are: 

Option  1:  Half  sine  wave  deceleration  pulse  Option  4:  Spline  fit  position,  velocity,  or 

Option  2:  Tabular  unidirectional  deceleration  deceleration  data. 

Option  3:  Six-degree-of-freedom  deceleration 


Option  1,  as  depicted  by  Figure  30,  is  the  simplest  way  to  define  vehicle  motion.  The  half  sine  wave 
in  the  top  graph  is  obtained  by  inputting  an  initial  velocity  and  time  duration,  VTIME,  of  the  pulse. 
The  magnitude  of  the  deceleration  is  automatically  calculated  so  that  the  final  velocity  is  zero.  The 
orientation  of  the  deceleration  vector  is  given  by  azimuth  and  elevation  angles.  A  set  of  X,  Y,  Z 
coordinates  is  used  to  specify  the  vehicle’s  initial  position  in  the  reference  system.  This  option  is  often 
used  to  approximate  the  vehicle  motion  when  detailed  motion  data  are  unavailable.  Example.ain  uses 
this  option.  Option  2  is  similar  to  option  1  except  that  the  half  sine  wave  is  replaced  by  a  table  of 
unidirectional  deceleration  data  supplied  at  a  fixed  time  interval,  ATD,  as  in  the  bottom  graph. 


Figure  30.  Prescribed  Motion  Options  1  and  2 
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To  define  more  complex  three-dimensional  vehicle  motion,  option  4  is  recommended  for  its  spline  fit 
capability.  Figure  3 1  demonstrates  the  input  involved  in  this  option.  A  table  of  data  of  the  vehicle’s 
position,  velocity,  or  acceleration  vs.  time  is  used  to  generate  spline  fit  functions.  If  the  table  is  a 
position  table,  the  position  vector  (X,^,  Xy,  XJ  and  the  vehicle’s  yaw,  pitch,  and  roll  (62,0y,  6J  are 
input.  If  it  is  a  velocity  table,  the  velocity  vector  (v^,  Vy,  v^)  and  angular  velocity  vector  (to^,  Wy,  (o^ 
are  input.  An  acceleration  table  uses  the  deceleration  vector  (a^,  ay,  a^)  and  angular  acceleration  vector 
(a^,  Oy ,  cQ.  The  ATB  Model  will  spline  fit  these  data  using  polynomials  and  then  compute  the  vehicle 
motion  at  a  fixed  time  interval,  set  by  the  user.  In  spline  fitting  angular  positions,  the  angle  values  will 
be  transformed  into  quaternions  and  the  four  quaternion  components  are  each  spline  fit  independently. 
The  yaw,  pitch,  and  roll  are  then  computed  from  the  quaternions.  Option  3  is  similar  to  option  4 
except  that  only  decelerations  are  allowed  and  the  data  must  be  supplied  at  even  time  points.  The 
decelerations  are  used  directly  by  the  program  and  spline  fitting  is  not  required. 


2.5.2  Initial  Positioning  of  Body 

For  body  segments,  the  initial  positions  and  velocities  are  defined  by  a  set  of  G  cards.  First,  a  G.2  card 
is  used  to  specify  the  initial  velocity  and  CG  position  of  each  body’s  reference  segment  in  a  vehicle 
or  groirnd  (inertial)  reference  system.  Initial  segment  velocities  can  be  set  equal  to  a  vehicle’s  initial 
velocity.  Once  the  reference  segment  is  positioned,  the  initial  yaw,  pitch,  roll,  and  angular  velocities 
of  all  the  body  segments’  local  reference  systems  are  specified  by  a  set  of  G.3  cards.  This  in  effect 
defines  all  the  segments’  initial  conditions  since  the  body  is  modeled  in  a  chain  structure.  Therefore, 
in  the  process  of  setting  up  a  body’s  initial  conditions,  the  positioning  of  the  reference  segment  in  the 
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G.2  cards  defines  the  whole  body’s  location  and  the  specifying  of  angles  and  angular  velocities  in  the 
G.3  cards  defines  the  body’s  posture.  These  initial  orientations  and  velocities  may  be  specified  with 
respect  to  the  ground,  any  vehicle,  or  any  segment  previously  positioned.  Figure  32  depicts  this 
concept. 


For  a  body  free  in  space,  the  above  process  can  be  relatively  straightforward;  however,  when  the  body 
is  confined  by  the  surrounding  environment  (for  example,  sitting  on  a  seat),  the  process  can  become 
fairly  involved.  The  reason  for  this  is  that  the  body  must  initially  be  in  static  equilibrium,  and  this 
equilibrium  is  achieved  by  balancing  gravitational  forces  against  contact  forces.  The  latter  are  highly 
position-dependent  and  must  be  properly  chosen  to  avoid  large  initial  segment  accelerations  due  to 
unbalanced  forces  and  torques  on  the  body.  An  iterative  adjustment  process  is  often  used  to  achieve 
this  objective.  Li  adjusting  the  segment  orientations,  care  should  be  taken  that  the  joints  are  not 
positioned  beyond  their  stops  or  inconsistent  with  their  rotational  constraints. 

The  iterative  method  requires  that  the  simulation  is  executed  to  zero  time.  A  tabular  printout  of  all  of 
the  external  forces  and  torques  and  resulting  linear  and  angular  accelerations  is  produced  for  time  zero. 
Joint  torques  can  also  be  checked  for  inconsistencies.  Then  the  user  adjusts  the  positions  based  on  the 
contact  forces  and  the  initial  linear  and  angular  accelerations.  This  procedure  usually  requires  several 
iterations  to  ensure  that  the  body  is  in  static  equilibrium  with  its  environment,  which  is  determined  by 
the  absence  of  large  accelerations  for  any  of  the  body  segments.  Perfect  equilibrium  is  generally  not 
attainable  for  the  seated  or  standing  position;  however,  small  initial  accelerations  are  tolerable, 
especially  for  the  angular  accelerations  which  are  usually  difficult  to  balance,  as  long  as  they  are  much 
smaller  than  the  accelerations  induced  by  the  dynamic  environmental  conditions  under  study. 
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For  initially  balancing  seated  human  or  dummy  subjects,  the  following  steps  are  suggested: 

1.  Start  position. 

-  Torso  segments  parallel  with  the  seat  back. 

-  Legs  rotated  above  the  floor  and  pedals. 

2.  Adjust  the  lower  torso  position  to  balance  more  than  60%  of  the  total  body  weight. 

-  Rotate  torso  segments  to  contact  seat  back. 

3.  Rotate  upper  legs  down  to  balance  most  of  the  body  weight. 

-  Check  horizontal  forces.  The  combined  forces  should  be  close  to  zero. 

4.  Rotate  lower  legs  to  fully  balance  weight  with  foot/floor  contacts. 

2.6  Dimension  Units,  Gravity,  and  Time  Control  of  the  ATB  Program 

Before  any  body  or  vehicle  data  can  be  considered  for  input  to  the  model,  the  user  must  decide  the 
units  of  measurement  to  be  used  and  the  gravity  direction,  relative  to  the  inertial  system. 

2.6.1  Selecting  Dimensional  Units 

The  units  of  measurement  for  the  input  data  (i.e.,  pounds/inches/seconds  or  Newtons/meters/seconds) 
must  be  chosen.  The  choice  is  arbitrary  and  there  is  no  default,  but  once  the  selection  is  made,  all 
input  data  must  be  in  the  same  units.  Choosing  the  units  of  measurement  for  the  input  data  also 
automatically  specifies  the  units  for  the  output  data.  The  units  of  measurement  are  selected  by 
supplying  the  alphanumeric  names  of  the  abbreviations  for  the  units  of  force  (UNITM),  distance 
(UNITL)  and  time  (UNITT). 

The  units  of  measurement  used  in  the  ATB  Input  Description  for  illustrative  purposes  are  pounds, 
inches,  and  seconds.  These  units  were  selected  at  the  time  of  the  initial  development  of  the  model 
when  most  available  data  were  in  these  units.  Although  there  are  no  official  units,  the  format  (field 
width  and  number  of  digits  following  the  decimal)  for  various  input  and  output  items  was  established 
on  the  basis  of  the  expected  magnitude  of  these  data  for  a  simple  car  crash  type  simulation,  assuming 
the  pound,  inch,  and  second  measurement  system.  Hence  it  is  possible  that  a  different  choice  of  units 
may  result  in  data  that,  while  numerically  correct,  may  not  fit  in  the  specified  input  and  output  format. 

Note  that  mass  units  are  not  required  for  input  and  output  purposes,  although  they  are  assumed 
internally  by  the  program.  This  is  accomplished  by  supplying  the  weight  of  the  body  segments  using 
the  force  units,  ibe  ATB  program  converts  these  input  values  to  mass  units  by  dividing  these  force 
units  by  the  value  of  the  acceleration  due  to  gravity,  which  must  be  provided  as  input.  Unfortunately, 
an  inconsistency  was  introduced  during  the  early  development  of  the  program  for  the  principal 
moments  of  inertia  input  units,  hi  retrospect,  the  units  for  the  principal  moments  of  inertia  should  have 
been  weight  (force)  multiplied  by  distance  squared,  and  the  input  values  converted  by  the  program  by 
dividing  by  the  acceleration  due  to  gravity,  as  is  done  for  the  segment  weights.  As  the  input  is  now 
established,  the  required  units  for  these  principal  moments  of  inertia  are  weight  (force)  multiplied  by 
distance  multiplied  by  time  squared,  which  is  equivalent  to  mass  multiplied  by  distance  squared.  This 
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inconsistency  has  never  been  removed  because  its  removal  would  invalidate  many  already  established 
input  files. 

Finally,  both  the  GEBOD  and  VIEW  programs  employ  similar  unit  conventions  as  the  ATB  Model. 

2.6.2  Specifying  Gravity 

Once  the  units  of  measurement  have  been  selected,  the  user  must  define  what  is  meant  by  the  inertial 
coordinate  system.  As  discussed  earlier,  the  inertial  coordinate  system  is  the  coordinate  system  to 
which  all  other  coordinate  systems  are  referred  and  it  is  within  this  system  that  Newton’s  laws  hold. 
The  inertial  coordinate  system  of  the  model  is  assumed  to  be  at  rest,  but  is  designated  as  a  segment 
called  the  ground  segment  with  its  segment  number  given  by  NGRND.  The  inertial  coordinate  system 
is  defined  by  specifying  the  gravity  vector  GRAVTY.  Most  whole-body  simulations  define  the  gravity 
vector  GRAVTY  as  (0, 0,  g),  to  be  aligned  with  the  positive  Z  axis  of  the  inertial  coordinate  system, 
g  corresponds  to  the  standard  coefficient  of  gravity  at  the  surface  of  the  earth. 

The  gravity  field  defined  by  GRAVTY  is  assumed  to  be  constant  throughout  space  and  time  in  the 
ATB  Model  and  is  applied  to  all  segments  that  are  given  a  nonzero  weight.  The  magnitude  of  the 
vector  GRAVTY  is  used  to  compute  the  masses  of  the  segments  from  their  supplied  weights.  If  the 
user  wants  to  simulate  the  motion  of  an  object  in  a  zero  gravity  field,  such  as  a  spacecraft  in  deep 
space,  the  gravity  vector  would  be  supplied  as  GRAVTY  (0,  0,  0).  The  magnitude  of  this  vector  is 
obviously  zero,  so  computation  of  the  masses  of  the  segments  from  their  weights  would  not  be  possible 
using  the  magnimde  of  GRAVTY.  To  circumvent  this  problem,  the  user  has  the  option  of  supplying 
G.  G  represents  a  factor  by  which  the  weights  of  the  segments  will  be  divided  to  yield  a  mass.  If  G 
is  supplied  as  nonzero,  the  ATB  program  will  use  the  value  of  G  (rather  than  the  magnitude  of 
GRAVTY)  to  compute  the  masses  of  the  segments  and  will  apply  the  force  vector,  GRAVTY,  to  all 
segments  with  a  nonzero  mass.  G  must  be  nonzero  when  GRAVTY  (0, 0, 0)  is  used. 

2.6.3  Integration  and  Output  Time  Control 

Time  control  parameters  must  be  specified  for  each  simulation.  These  parameters  control  the  length 
(simulation  time)  of  the  run,  the  amount  and  format  of  the  output,  the  tabular  time  histories,  and 
operation  of  the  program  integrator.  Although  the  program  places  no  restrictions  on  these  input 
parameters,  a  judicious  choice  of  the  parameters  can  improve  computational  efficiency  and  numerical 
stability.  For  the  purpose  of  clarity,  the  following  description  will  assume  that  UNITT,  units  of  time, 
are  seconds. 

Figure  33  illustrates  the  time  control  mechanism.  After  all  input  and  initialization  are  performed,  time 
is  advanced  in  overall  steps  of  DT  seconds.  The  optional  output  is  provided  at  time  zero  and  each 
integral  multiple  of  DT  seconds  of  simulation  time.  The  total  simulation  time  is  NSTEPS*DT 
seconds  where  NSTEPS  and  DT  are  input  parameters.  The  values  of  NSTEPS  and  DT  should  be 
chosen  to  provide  the  desired  length  of  the  simulation,  and  the  amount  and  frequency  of  output  data. 
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A  secondary  control  of  time  is  performed  by  the  program  integrator.  Integrator  steps  are  controlled 
by  the  supplied  values  for  HO,  HMIN,  and  HMAX.  The  integrator  advances  time  in  substeps  of  H 
seconds  starting  with  HO  and  varying  between  HMIN  and  HMAX.  H  is  halved  when  convergence 
of  the  integrator  parameters  is  not  achieved,  but  H  is  not  permitted  to  become  less  than  HMIN.  K 
convergence  is  not  attained  with  an  HMIN  time  step,  the  simulation  stops  with  a  descriptive  error 
message  written  to  the  .aou  file.  If  the  convergence  criteria  are  satisfied  for  several  integration  steps, 
H  is  doubled.  This  can  continue  until  the  integration  step  reaches  HMAX.  H  is  also  adjusted  to 
provide  an  integrator  time  step  at  each  integer  multiple  of  DT,  so  that  output  data  is  available.  For  this 
reason,  it  is  best  that  DT  be  an  integer  multiple  of  HMAX  so  that  the  DT  time  step  will  be  executed 
in  equal  HMAX  substeps  during  periods  of  stable  activity.  Also,  since  the  value  of  H  is  permitted  to 
double  during  these  stable  periods  or  be  halved  during  unstable  periods,  the  integrator  will  execute 
more  efficiently  if  HMAX  is  a  power  of  two  multiple  of  HMIN  and  HO. 


tsO  — - - ■>  t=DTxNSTEP 


It  has  been  observed  that  suitable  values  for  HMAX  lie  between  one  and  five  msec  for  most  occupant 
and  pedestrian  simulations.  Generally,  values  for  DT  of  0.002, 0.004, 0.010,  or  0.020  seconds;  for 
HMAX  of  0.001  or  0.002  seconds;  and  for  HMIN  and  HO  of  0.000125  or  0.000250  seconds  work 
satisfactorily.  It  is  possible  to  execute  the  integrator  in  a  “fixed  step  mode”  by  setting  HMAX  = 
HMIN  =  HO,  but  this  is  not  recommended. 
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3.  ORGANIZATION  OF  ATB  INPUT  DATA  AND  CONTROL  OF  OUTPUT  DATA 


3.1  Structure  of  Input  File 

The  input  for  the  ATB  program  is  contained  in  a  single  primary  input  file  (FORTRAN  unit  No.  5). 
It  is  a  formatted  file,  structured  in  a  fixed  80  column  card  format,  of  alphanumeric  data  input  (see 
example.airi).  Each  record  or  line  of  the  file  therefore  corresponds  to  the  contents  of  an  input  card  that 
has  a  unique  identification  (e.g.,  input  Card  A.l.a).  This  produces  a  modular  form  for  the  contents  of 
an  input  file  for  the  ATB  program.  For  example,  the  A  input  cards  contain  the  general  ran  parameters, 
the  B  input  cards  contain  the  inertial  and  geometric  parameters  that  define  the  segments  and  joints  of 
the  body,  the  C  input  cards  contain  the  parameters  that  define  the  prescribed  motions,  etc. 

During  the  input  portion  of  the  ATB  program  execution,  considerable  program  initialization  is 
performed  and  a  completely  annotated  listing  of  the  program  input  is  produced  on  the  primary  output 
unit  (FORTRAN  unit  No.  6),  as  shown  at  the  beginning  of  example.aou. 

Following  is  a  summary  of  all  of  the  input  cards.  A  complete  description  giving  the  format  for  each 
card,  the  conditions  that  specify  its  necessity,  the  input  parameters  to  be  supplied  on  each  card,  and 
a  definition  of  each  of  these  parameters  is  provided  in  the  ATB  Model  Input  Manual. 

Card  A.  Run  control  parameters 

A.  1  .a-c  Date  and  comment 

A.2  Not  used 

A.3  Dimensional  units,  components  of  gravity 

A.4  Integrator  parameters 

A. 5  NPRT  array  for  output  control 

Card  B.  Physical  characteristics  of  the  body 

B.  1  Body  title,  number  of  segments,  joints,  and  of  deformable  segments 

B.  1  .b-e  Finite  element  analysis  data  for  deformable  segments 

B.2.a-b  Physical  characteristics  of  body  segments 

B.3.a-b  Physical  characteristics  of  joints 

B.3.C  Node  numbers  used  to  compute  rotational  deformation  for 

deformable  segments 

B.4  Joint  spring  function  coefficients,  or  restoring  torque  function  numbers 

B.5  Joint  viscous  function  coefficients 

B.6  Integrator  convergence  tests  for  body  segments 

B.7.a-b  Controls  for  flexible  elements 


38 


CardC.  Prescribed  motion 


C.  1  Vehicle  motion  title 

C.2.a-b  Prescribed  motion  control  parameters 

C.3  Unidirectional  deceleration  tables 

C.4  Six  degrees-of-freedom  deceleration  tables 

C. 5  Spline  fit  tables 

Card  D.  Contact  surface  and  other  environment  definitions 

D.  1  .a-b  Number  of  contact  panels,  belts,  airbags,  etc.,  and  water  force  and  joint 

actuator  switches 

D.2.a-d  Plane  description  and  input  data 

D.3.a-c  Simple  belt  description  and  input  data 

D.4.a-h  Airbag  description  and  input  data 

D.5  Additional  contact  (hyper)ellipsoid  data 

D.6  Constraint  and  tension  element  input  data 

D.7  Body  segment  symmetry  options 

D.  8 .  a  Spring-damper  input  data 

D.8.b  Spring-damper  attachment  node  for  deformable  segments 
D.9.a  Applied  force/torque  function  input  data 

D. 9.b  Deformable  segments’  node  numbers  for  applied  force/torque 

Card  E.  Function  definitions 

E.  1  Function  identification  number  and  title 

E.2  Function  definition  control  parameters 

E.3  5th  degree  polynomial  coefficients 

E.4.a-b  Tabular  function  definition 

E.5  (No  longer  required  by  program) 

E.6.a-d  Wind  force  functions  input  data 

E. 7.a-d  Joint  restoring  force  functions  input  data 

Card  F.  Allowed  contacts  and  associated  functions 

F.  1  .a-b  Plane/ellipsoid  contact  definition 

F.2.a-b  Belt/segment  contact  definition 

F.3.a-b  Ellipsoid/ellipsoid  contact  definition 

F.4.a-b  Specifications  for  globalgraphic  joint  functions 

F.5.a  Not  used 

F.6  Airbag/segment  contact  definition 

F.7.a-b  Wind  force  function  specification 

F.8.a-d  Hamess/belt  system  input  data 


F.9.a-m 

F.IO 


Water  force  simulation  data  input 
Joint  actuator  data  input 


Card  G.  Initial  positioning  input 


G.l.a 

Segment  initial  velocity  data  source 

G.l.b 

Not  used 

G.2 

Initial  position  and  velocity  for  reference  segments 

G.3.a-b 

Initial  segment  angular  orientation  and  velocity  input  data 

G.4 

Equilibrium  control  parameters 

G.5 

Equilibrium  control  assignments 

G.6 

Equilibrium  constraint  assignments 

Card  H.  Tabular  time  history  output  control  parameters 

H.l.a-b 

Linear  accelerations  of  selected  points  on  segments 

H.2.a-b 

Linear  velocities  of  selected  points  on  segments 

H.3.a-b 

Linear  positions  of  selected  points  on  segments 

H.4 

Angular  accelerations  of  selected  segments 

H.5 

Angular  velocities  of  selected  segments 

H.6 

Angular  orientations  of  selected  segments 

H.7 

Joint  parameters  for  selected  joints 

H.8 

Wind  forces  on  selected  segments 

H.9 

Joint  forces  and  torques  for  selected  joints 

H.10.a-c 

Properties  of  selected  groups  of  segments 

H.ll 

Actuator  joint  torques 

H.12 

Parameters  for  HIC,  HSI,  and  CSI  computations 

3.2  ATB  Output  Files 

Because  of  the  complexity  of  the  ATB  Model  and  the  potential  for  huge  amounts  of  output  from  a 
single  simulation,  the  ATB  program  was  written  so  that  possible  output  files  are  controlled  for  each 
run.  This  tailoring  of  the  number  of  files  to  be  written  for  each  simulation  has  been  somewhat 
confusing  because  some  aspects  of  it  are  explicit  (the  user  sets  a  flag  for  the  type  and  frequency  of  the 
desired  output)  and  others  are  implicit  (indirectly  determined  by  the  type  and  number  of  force 
deflection  interactions,  etc). 

A  logical  unit  is  the  device  or  file  from  which  or  to  which  input  or  output  from  a  FORTRAN  program 
is  to  be  sent.  Except  for  the  primary  input  and  output  files  (FORTRAN  unit  Nos.  5  and  6),  the  use  of 
each  I/O  file  is  controlled  by  input  parameters  contained  within  the  program  input  file.  The  ATB 
model  has  an  open-ended  number  of  required  logical  units  which  depends  on  the  amount  of  output 
requested  by  parameters  in  the  input  file.  Table  2  summarizes  major  FORTRAN  logical  units  that 
may  be  used  by  the  ATB  program. 
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Table  2  Summary  of  ATB  Program  I/O  Units 


Logical 

Unit 

Filename 

Extension 

Type* 

Description 

Controlling  Parameters 

1 

.SAl,  .TPl,  .UFl 

UorF 

Program  VIEW  input 

NPRT(1)&(35)  on  CardA.5 

5 

.AIN 

F 

Primary  input 

always  required 

6 

.AOU 

F 

Primary  output 

always  required 

21+ 

T?? 

F 

Time  histories 

NPRT  (4)  on  Card  A.5  &  Cards  H.l~  H.12 

*  Type  is  F  for  formatted,  U  for  unformatted  file. 


3.2.1  View  Output  (Unit  1) 

Logical  unit  No.  1  is  typically  an  ASCII-formatted  output  file  designed  to  be  used  as  data  input  to 
VIEW,  the  program  that  creates  the  graphics  frequently  associated  with  the  ATB  Model.  This  output 
file  has  an  extension  name  .sal,  oi.tpl.  For  backward  compatibility  with  an  early  version  of  VIEW, 
an  optional  unformatted  (binary)  output  file  is  also  included,  with  an  extension  name  .ufl. 

The  generation  of  this  output  file  is  controlled  by  the  value  of  NPRT  (1)  that  is  supplied  on  input  Card 
A.5.  A  blank  or  zero  value  for  NPRT  (1)  will  suppress  the  generation  of  output  file  .sal  or  .tpl, 
whereas  a  non-zero  positive  value  will  produce  data  records  that  are  equally  spaced  at  every  m*DT 
seconds  of  simulation  time  starting  at  0  time,  where  m  is  the  integer  value  of  NPRT  (1)  and  DT  is 
defined  on  input  Card  A.4. 

File  .sal  and  .tpl  contain  fixed  initialization  data  describing  the  planes,  contact  ellipsoids,  and  harness 
belts.  These  are  followed  by  records  containing  the  values  of  time  and  the  corresponding  dynamic 
data,  including  segment  positions  in  the  inertial  reference  and  the  direction  cosine  matrix  for  each  of 
the  body  segments.  The  only  difference  between  the  .sal  and  .tpl  files  is  the  formatting.  The  .sal  file 
uses  the  new  structured  ASCII  graphics  output  format  used  by  VIEW  and  is  selected  by  setting 
NPRT(35)  to  0.  The  structured  graphics  output  file  is  designed  for  easier  troubleshooting  and  smaller 
file  size.  The  .tpl  file  uses  the  old  ASCII  graphics  output  format  used  by  older  versions  of  VIEW  and 
is  selected  by  setting  NPRT(35)  to  2. 


41 


3.2.2  Primary  Output  (Unit  6) 

The  primary  output  file  for  the  ATB  program  is  logical  unit  No.  6.  It  has  an  extension  name  .aou. 
Except  for  injury  criteria  (HIC,  HSI,  etc.)  results,  it  is  recommended  that  this  file  be  used  mainly  for 
diagnostics  and  input  reference  instead  of  simulation  result  analysis.  Referring  to  example,  aou  in  the 
Appendix,  the  primary  output  file  contains  the  following  items: 

1 .  A  labeled  echo  of  the  ATB  program  input  data. 

2.  Tables  of  segment  linear  and  angular  position,  velocity,  and  acceleration  information,  joint 
forces  and  torques,  the  sum  of  all  external  forces  and  torques  acting  on  each  segment,  and 
constraint  forces  data.  These  data  are  generated  at  fixed  time  intervals  of  m*DT  seconds, 
where  m  is  the  integer  value  of  NPRT(3).  These  tables  are  useful  in  determining  whether 
the  occupant  is  initially  balanced. 

3 .  Tables  of  the  computer  elapsed  CPU  time  used  by  selected  subroutines  and  the  number  of 
calls  to  these  subroutines.  They  are  printed  at  fixed  time  intervals  as  specified  by  DT  and 
at  a  frequency  specified  by  NPRT  (2).  When  NPRT  (2)  is  zero,  the  table  is  generated  only 
once  at  the  successful  completion  of  an  ATB  program  run. 

4.  Diagnostic-type  output  produced  at  every  call  to  various  subroutines,  as  controlled  by  the 
values  supplied  for  NPRT  (8)  to  (28)  on  input  Card  A.5.  This  output  is  intended  for 
diagnostic  or  checkout  purposes  only,  and,  if  used  indiscriminately,  can  produce  voluminous 
amounts  of  output. 

5.  Short  descriptions  of  changes  in  some  of  the  simulation  mn  conditions  are  produced  as  they 
occur.  They  include: 

a.  Failures  of  the  program  integrator  convergence  tests  that  cause  the  integration  step  to 
decrease  in  size.  The  time,  step  size,  segment  and  test  involved,  and  the  final 
convergence  test  parameters  are  printed.  NOTE:  These  messages  are  normal  and  do 
not  indicate  an  error  in  the  simulation.  A  stop  will  occur  if  the  integration  step 
becomes  too  small. 

b.  Changes  in  the  lock  conditions  of  joints  as  detected  by  changes  in  the  values  of  IPIN 
or  lEULER  for  the  various  joints.  The  time,  previous  and  new  values  of  the  indicator, 
and  the  identification  number  and  nomenclature  of  the  joint  involved  are  printed. 

c.  Each  time  a  point  is  added  to  or  deleted  from  the  set  of  harness  belt  reference  points, 
changes  are  indicated  by  listing  the  time,  the  set  of  points,  and  the  distance  between 
them. 
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6. 


A  page  containing  values  of  the  head  injury  criterion  (HIC),  head  and  chest  severity  indices 
(HSI  and  CSI),  and  related  information. 

7.  The  tabular  time  histories  may  be  generated  on  the  primary  output  file  as  described  in  the 
next  section. 

3.2.3  Tabular  Time  Histories  (Units  21, 22, 23, ...) 

The  tabular  time  histories  are  perhaps  the  most  useful  output  of  the  ATB  program.  Most  ATB 
simulation  result  analyses  are  performed  with  these  time  history  files.  These  files  are  designed  in  such 
a  way  as  to  allow  easy  data  export  to  spreadsheet  software,  such  as  Microsoft  Excel  and  SigmaPlot. 
Depending  on  the  value  of  NPRT(4),  they  may  be  output  at  the  end  of  the  .aou  file  or  to  consecutive 
logical  units  starting  from  UNIT21,  with  one  logical  unit  for  each  time  history.  For  the  latter,  each 
time  history  file  has  an  extension  name  .txx,  where  xx  is  the  corresponding  logical  unit  number.  Their 
generation,  contents,  frequency  of  output,  and  the  manner  by  which  they  are  generated  are  completely 
controlled  by  program  input  parameters.  The  Appendix  includes  several  example  time  history  files. 

There  are  two  types  of  time  history  files.  The  first  type  is  controlled  by  Cards  H.  1  to  H.  1 1 .  The  files 
in  this  set  primarily  contain  the  body  kinematic  data  and  properties.  The  output  of  each  individual  file 
in  this  set  is  optional  and  is  defined  by  the  user’s  specification  in  the  H  cards.  If  there  are  any  output 
files  of  this  type,  they  will  be  assigned  logical  unit  numbers  starting  with  21.  The  second  type  of 
output  files  include  contact  results,  such  as  deflections,  forces,  and  contact  locations,  and  spring- 
damper  forces.  Their  output  is  controlled  by  the  value  of  NPRT(18)  on  input  Card  A.5  and  an  output 
flag  in  each  contact  definition.  For  plane/ellipsoid  and  ellipsoid/ellipsoid  contacts,  the  value  of 
NPRT(18)  determines  whether  any  data  will  be  output.  If  NPRT(18)  allows  this  output,  then  the  user 
can  specify  which  individual  contacts  to  output  on  the  F.  1  and  F.3  cards.  This  gives  users  the  ability 
to  output  only  those  contacts  of  interest. 

Table  3  summarizes  the  time  history  files  and  their  control  cards.  The  listing  order  of  the  time 
histories  reflects  the  order  in  which  they  are  assigned  logical  unit  numbers. 
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Tables  Time  History  Files 


Time  History 

Output  Control  Cards 

Point  Linear  Accelerations:  Components  and  Resultant 

H.l 

Point  Linear  Velocities:  Components  and  Resultant 

H.2 

Point  Linear  Positions:  Components  and  Resultant 

H.3 

Segment  Angular  Accelerations:  Components  and  Resultant 

H.4 

Segment  Angular  Velocities:  Components  and  Resultant 

H.5 

Segment  Angular  Positions:  Yaw,  Pitch,  Roll,  and  Resultant 

H.6 

Joint  Parameters:  Lock  Condition,  Angles,  and  Resistive  Torques 

H.7 

Segment  Wind  Forces:  Components  and  Resultant 

H.8 

Joint  Forces  and  Torques:  Components  and  Resultant 

H.9 

Total  Body  Properties:  CG  Location,  Linear  and  Angular 
Momentum,  Kinetic  Energy,  Inertial  Tensor  Matrix,  Principal 
Moments  of  Inertia  and  Principal  Axes 

H.IO 

Active  (Actuator)  Joint  Torques:  Components  and  Resultant 

H.ll 

Plane/Ellipsoid  Contacts:  Normal,  Friction,  and  Resultant  Contact 
Forces,  Deflection,  and  Contact  Point  Coordinates 

A.5  andF.l 

Simple  Belt  Contacts:  Belt  Strain  and  Anchor  Point  Forces 

A.5 

Harness  Belt  Contacts:  Belt  Strain  and  Anchor  Point  Forces 

A.5 

Spring  Damper  Forces:  Component  and  Resultant 

A.5 

Ellipsoid/Ellipsoid  Contacts:  Normal,  Friction,  and  Resultant 
Contact  Forces,  Deflection,  and  Contact  Point  Coordinates 

A.5  andF.S 

Airbag  Contacts:  Airbag  Parameters  and  Contact  Forces 

A.5  and  F.6 
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APPENDIX  A 


ATB  Simulation  Example 

A  standard  sled  test  simulation  is  used  as  the  example  throughout  this  User’s  Guide.  In  this  example, 
a  167  lb  male  subject  was  the  occupant  restrained  using  a  double  shoulder  harness  and  a  lap  belt  with 
a  negative  G  strap.  A  generic  seat  was  used  with  the  seat  back  reclined  13  °  from  vertical  and  the  seat 
pan  inclined  6®  from  horizontal.  The  acceleration  waveform  of  the  impact  sled  was  an  approximate 
half-sine  pulse  with  a  209  msec  impact  duration  and  a  9.47  G  amplitude.  The  ATB  simulation  time 
was  300  msec.  Figure  A-1  shows  the  simulation  graphics  generated  by  the  VIEW  program.  VIEW 
reads  in  data  from  the  ATB  output  file  example.sal  and  draws  the  graphics.  Example.sal  is  not 
included  in  this  appendix  since  it  is  merely  a  data  file.  Corresponding  frames  from  high  speed  film 
are  presented  to  offer  a  comparison  between  the  ATB  simulation  results  and  the  actual  occupant 
response.  A  set  of  simulation  input  and  output  files  are  included  in  the  following  sections.  It  is 
recommended  that  this  example  be  used  as  a  trial  run  for  users  learning  the  ATB  program. 
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Figure  A-1.  Sled  Test  and  ATB  Simulation 
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Example.ain  File 


Example.ain  is  the  ATB  Model  input  file  for  the  sled  test  simulation.  Each  line  is  labeled  with  a  card 
number  at  the  end  of  the  line.  (The  term  “card”  is  a  carryover  from  the  days  of  punchcards.)  The  B 
Cards  are  generated  using  the  GEBOD  program.  The  user  should  refer  to  the  ATB  Model  Input 
Manual  for  detailed  variables  and  input  format  descriptions. 
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21  FEB  1997  0  00.000000  CARD  A1 

SIMULATION  OF  THE  HUMAN  VOLUNTEER  SLED  TEST  CARD  A2 

USING  NEWLY  DEVELOPED  HUMAN  JOINT  CHARACTERISTICS  CARD  A2 

IN.  LB. SEC.  0.000000  0.000000  386.0880  0.000000  CARD  A3 

4  150.0020000.0005000.0010000.0000625  CARD  A4 

555000000000000000100000  0-1  000000000  ICARD  A5 

1  5  1  15  CARD  A6 

15  14  MALE  HUMAN  167  LB  CARD  Bl 

LT  123. 597. 84590. 77400. 977304. 61706. 85004. 2530-. 4620. 00000. 83700  1  CARD  B2 

.00000.00000.00000 

CT  24.8720.13360.07370.200904.29105.97103.3030-1.430.00000.09500  1  CARD  B2 


.00000.00000.00000 

UT  350. 5914. 07732. 99622. 52424. 73806. 35807. 0970. 00000. 00000-. 1080  1  CARD  B2 

.0000014.400.00000 

N  42. 1790. 01450. 01750. 021602. 37602. 37604. 3700-. 4750. 000001. 1710  1  CARD  B2 

.00000.00000.00000 

H  59.2360.17970.20490.132803.90103.06105.6610-1.115.00000.00000  1  CARD  B2 

.0000036.000.00000 

RUL  620. 313. 25680. 25550. 014403. 00103. 001011. 623. 00000-. 3210. 50900  1  CARD  B2 

.00000.00000.00000 

RLL  78.0080.49270.50010.056902.29702.297010.188.91900-1.110.87200  1  CARD  B2 

.00000.00000.00000 

RF  82. 0100. 03840. 03640. 006901. 38301. 90905. 2870-.0230-. 6270. 00000  1  CARD  B2 

-4.0008.4000-6.100 

LUL  920.313.25680.25550.014403.00103.001011.623.00000.32100.50900  1  CARD  B2 

.00000.00000.00000 

T.T.T.  a8. 0080. 49270. 50010. 056902. 29702. 297010. 188. 919001. 1100. 87200  1  CARD  B2 

.00000.00000.00000 

LF  b2. 0100. 03840. 03640. 006901. 38301. 90905. 2870-. 0230. 62700. 00000  1  CARD  B2 

4.00008.40006.1000 

RUA  C4. 0180. 10680. 11310. 021601. 93901. 93906. 5350. 00000-.2230-. 5490  1  CARD  B2 

.00000.00000.00000 

RLA  D3. 9140. 25680. 25550. 014401. 74201. 74209. 1950. 00000. 611001. 1510  1  CARD  B2 

.00000.00000.00000 

LUA  E4. 0180. 10680. 11310. 021601. 93901. 93906. 5350. 00000. 22300-. 5490  1  CARD  B2 

.00000.00000.00000 

LLA  F3. 9140. 25680. 25550. 014401. 74201. 74209. 1950. 00000-. 61101. 1510  1  CARD  B2 

.00000.00000.00000 

P  1  0-1.400.00000-2.230-2.340.000002.2100  0.00000.00000  CARD  B3 

.00000.00000.00000.000005.0000.00000.00000.00000.00000  321321 
W  2  0-1. 680. 00000-.8300-. 2700. 000006. 9900  0.00000.00000  CARD  B3 

.00000.00000.00000.000005.0000.00000.00000.00000.00000  321321 
NP  3  0-. 2200. 00000-7. 210-. 8800. 000001. 4900  0.00000.00000  CARD  B3 

.00000.00000.00000.0000010.000.00000.00000.00000.00000  321321 
HP  4  0.98000. 00000-2. 500-. 8200. 000001. 9600  0.00000.00000  CARD  B3 

.00000.00000.00000.0000010.000.00000.00000.00000.00000  321321 
RH  1  0-. 54002. 11001. 4800-. 4100-1. 920-7. 500  0.00000.00000  CARD  B3 

14.00048.000.00000.00000.00000-7.000.00000.00000.00000  321321 
PUC  6  1-. 2400. 360009. 7000. 63000-. 5500-6. 890  0.00000.00000  CARD  B3 

.00000-66.00.00000.00000.00000.00000.00000.00000.00000  321321 
RA  7  0.39000-. 75009. 67001. 3600-. 3100-2. 740  0.00000.00000  CARD  B3 

.0000065.000.00000.00000.00000.00000.00000.00000.00000  321321 
LH  1  0-. 5400-2. 1101. 4800-. 41001. 9200-7. 500  0.00000.00000  CARD  B3 

-14.0048.000.00000.00000.000007.0000.00000.00000.00000  321321 
LK  9  1-.2400-. 36009. 7000. 63000. 55000-6. 890  0.00000.00000  CARD  B3 

.00000-66.00.00000.00000.00000.00000.00000.00000.00000  321321 
LA  10  0.39000.750009.67001.3600.31000-2.740  0.00000.00000  CARD  B3 

.0000065.000.00000.00000.00000.00000.00000.00000.00000  321321 
RS  3  0-. 93006. 5200-4. 240. 51000-. 2200-5. 150  0.00000.00000  CARD  B3 

59.29079.080.00000.00000.00000.00000.00000.00000.00000  321321 
RE  12  1-.6500-. 39005. 0200-. 5000. 31000-7. 080  0.00000.00000  CARD  B3 

-15.0065.000.00000-15.00.00000.00000.00000.00000.00000  321321 
LS  3  0-. 9300-6. 520-4. 240. 51000. 22000-5. 150  0.00000.00000  CARD  B3 

-59.2979.080.00000.00000.00000.00000.00000.00000.00000  321321 
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1-. 6500. 390005. 0200-.5000-. 3100-7. 080  0.00000.00000  CARD 

15.0065.000.00000  15.00.00000.00000.00000.00000.00000  3213 


.0000010.000.00000.70000 

.0000010.000.00000.70000 

.000004.0000.00000.70000 

.000004.0000.00000.70000 


-4  .00000.00000.00000  0. 
-10  1.8000.00000.70000  60 
-12  7.0000.00000.70000  35 
-5  .00000.00000.00000  0. 
-11  1.8000.00000.70000  60 
-15  7.0000.00000.70000  35 
-1  .00000.00000.00000  0. 
-8  1.8000.00000.70000  70 
-2  .00000.00000.00000  0. 
-9  1.8000.00000.70000  70 
.50000.0000030.000.00000.00000 
.50000.0000030.000.00000.00000 


20.00000.0000010.000.00000.70000 
20.00000.0000010.000.00000.70000 
35.00000.0000010.000.00000.70000 
25 . 00000 . 0000010 . 000 . 00000 .70000 
0.000000. 00000. 20020-8E-04. 70000 
60.00000.00000.00000.00000.00000 
35.000005.8378.05646.00000.70000 
0.000000. 00000. 20020-8E-04. 70000 
60.00000.00000.00000.00000.00000 
35.000005.8378.05646.00000.70000 
0 . 000000 . 00000 . 07450-3E-04 .70000 
70.00000.00000.00000.00000.00000 
0.000000. 00000. 07450-3E-04. 70000 
70.00000.00000.00000.00000.00000 
00  .00000.00000 
00  .00000.00000 


5.000000CARD 

35.00000CARD 

35.00000CARD 

35.00000CARD 

lO.OOOOOCARD 

O.OOOOOOCARD 

4.000000CARD 

lO.OOOOOCARD 

O.OOOOOOCARD 

4.000000CARD 

25.00000CARD 

O.OOOOOOCARD 

25.00000CARD 

O.OOOOOOCARD 

CARD 


.60000.0000030.000.00000.00000 

.60000.0000030.000.00000.00000 

.30000.0000030.000.00000.00000 

.30000.0000030.000.00000.00000 

.30000.0000030.000.00000.00000 

.30000.0000030.000.00000.00000 

.30000.0000030.000.00000.00000 

.30000.0000030.000.00000.00000 

.60000.0000030.000.00000.00000 

.30000.0000030.000.00000.00000 

.60000.0000030.000.00000.00000 

.30000.0000030.000.00000.00000 


.00000.00000 

.00000.00000 

.00000.00000 

.00000.00000 

.00000.00000 

.00000.00000 

.00000.00000 

.00000.00000 

.00000.00000 

.00000.00000 

.00000.00000 

.00000.00000 


.01000. 01000. 01000. 01000. 01000. 01000. 10000. 10000. 10000. 10000. 10000. OlOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
.01000. 01000. 01000. 00000. 00000. 00000. 10000. 10000. 10000. 00000. 00000. OOOOOCARD 
SLED  ACCELERATION  CARD 

.00000.00000.00000.00000.00000.00000.00000.00000  57.00000.00400  OCARD 

.02560. 09520. 00140. 02560. 00140. 02560. 09520. 00140. 02560. 00140. 02560. 04690CARD 
.02270. 11930. 21600. 554201. 76202. 22102. 75243. 67044. 61255. 65126. 59337. 4629CARD 
8. 38098. 98489. 39559. 39559. 41969. 25059. 17819. 17818. 98488. 81578. 62258. 6225CARD 
8. 57418. 40508. 13937. 77707. 41467. 00406. 56916. 11025. 72375. 14394. 46753. 8878CARD 
3.25972.55921.8103.98900.26430.55700.50870.24300.00000  CARD 

12  0  0  0  0  0.1  0  0  0  CARD 

1  SEAT  CUSHION  0  0  CARD 


10.00000 

8.000000 

-10.00000 

CARD  D2B 

28.01000 

8.000000 

-11.89000 

CARD  D2C 

10.00000 

-8.000000 

-10.00000 

CARD  D2D 

SEAT  BACK 

0  0 

CARD  D2A 

1.000000 

9.000000 

-48.97000 

CARD  D2B 

10.00000 

9.000000 

-10.00000 

CARD  D2C 

1.000000 

-9.000000 

-48.97000 

CARD  D2D 

FLOOR. 

0  0 

CARD  D2A 

0.000000 

12.00000 

-1.300000 

CARD  D2B 

60.00000 

12.00000 

-1.300000 

CARD  D2C 
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0.000000 

-12.00000 

-1.300000 

4  HEAD 

PAD 

0  0 

2.891687 

7.499822 

-45.99358 

5.371688 

7.499822 

-35.28358 

2.891687 

-7.500182 

-45.99358 

5  SEAT 

FRONT  PANEL. 

0  0 

28.01000 

8.000000 

-11.89000 

26.66000 

8.000000 

-4.400000 

28.01000 

-8.000000 

-11.89000 

6  BACK 

PANEL2.  13  DEGR  0  0 

1.000000 

9.000000 

-48.97000 

10.00000 

9.000000 

-10.00000 

1.000000 

-9.000000 

-48.97000 

7  FIREWALL. 

0  0 

60.00000 

12.00000 

-25.00000 

60.00000 

-12.00000 

-25.00000 

60.00000 

12.00000  - 

•0.7500000 

8  RIGHT 

‘  SIDE  SEAT/ IN. 

0  0 

8.410000 

8.100000 

-6.660000 

8.700000 

8.100000 

-14.73000 

30.58001 

.  8.100000 

-6.640000 

9  LEFT 

SIDE  SEAT/ IN. 

0  0 

8.410000 

-8.100000 

-6.660000 

30.58001 

-8.100000 

-6.640000 

8.700000 

-8.100000 

-14.73000 

10  RUDDER  PEDALS. 

0  0 

44.99118 

8.999872 

-1.423467 

48.27194 

8.999872 

-3.562127 

44.99118 

-9.000128 

-1.423467 

11  LEFT 

SIDE  PANEL. 

0  0 

1.000000 

-9.000000 

-48.97000 

10.90000 

-9.000000 

-6.100000 

-7.770000 

-9.000000 

-46.95000 

12  RIGHT 

'  SIDE  PANEL. 

0  0 

1.000000 

9.000000 

-48.97000 

-7.770000 

9.000000 

-46.95000 

10.90000 

9.000000 

-6.100000 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

3  SEGMENT- SEGMENT  FCN.  0  0 

0.000000 

-5.000000 

0.000000 

0.000000 

6 

0.000000 

0.000000 

1.000000 

470.0000 

3.000000 

1220.000 

4.000000 

1470.000 

7  R  FACTOR. 

0  0 

0.000000 

0.000000 

0.7000000 

0.000000 

13  STIFF 

•  SURFACES 

0  0 

0.000000 

-4.000000 

0.000000 

0.000000 

8 

0.000000 

0.000000 

0.1000000 

5.000000 

0.3000000 

40.00000 

0.4000000 

60.00000 

2.000000 

1200.000 

3.000000 

4000.000 

14  FRICTION  FUNC . 

0  0 

0.000000 

0.000000 

2.0000000 

0.000000 

19  CF=.25,CREST=.25 

0  0 

0.000000 

0.000000 

0.2500000 

0.000000 

20  DAMPING  COEFF.  C=100  0  0 

0.000000 

1,000000 

0.000000 

0.000000 

0.000000 

1000.000 

0.000000 

0.000000 

21  RATE 

OF  DEFLEC. 

0  0 

-40.00000 

-150.0000 

0.000000 

0.000000 

21 

-40.00000 

0.000000 

-30.00000 

0.000000 

-10.00000 

0.000000 

0.000000 

0.000000 

10.00000 

1.000000 

20.00000 

0.9899999 

0  0  0 

1.000000 

2.000000 

5.000000 

0.000000 

1.000000 

0.2000000 

1.000000 


2.000000 

0.000000 

1.000000 

0.000000 

1.000000 

-20.00000 

5.000000 

30.00000 


CARD 

D2D 

CARD 

D2A 

CARD 

D2B 

CARD 

D2C 

CARD 

D2D 

CARD 

D2A 

CARD 

D2B 

CARD 

D2C 

CARD 

D2D 

CARD 

D2A 

CARD 

D2B 

CARD 

D2C 

CARD 

D2D 

CARD 

D2A 

CARD 

D2B 

CARD 

D2C 

CARD 

D2D 

CARD 

D2A 

CARD 

D2B 

CARD 

D2C 

CARD 

D2D 

CARD 

D2A 

CARD 

D2B 

CARD 

D2C 

CARD 

D2D 

CARD 

D2A 

CARD 

D2B 

CARD 

D2C 

CARD 

D2D 

CARD 

D2A 

CARD 

D2B 

CARD 

D2C 

CARD 

D2D 

CARD 

D2A 

CARD 

D2B 

CARD 

D2C 

CARD 

D2D 

CARD 

D7 

CARD 

El 

CARD 

E2 

CARD 

E4A 

889.9999CARD 

E4B 

1580.000CARD 

E4C 

CARD 

El 

CARD 

E2 

CARD 

El 

CARD 

E2 

CARD 

E4A 

20.00000CARD 

E4B 

560.0001CARD 

E4C 

CARD 

E4D 

CARD 

El 

CARD 

E2 

CARD 

El 

CARD 

E2 

CARD 

El 

CARD 

E2 

O.OOOOOOCARD 

E3 

CARD 

El 

CARD 

E2 

CARD 

E4A 

O.OOOOOOCARD 

E4B 

l.OOOOOOCARD 

E4C 

0.9650000CARD 

E4D 
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40.00000 

0.9279998 

50.00000 

0.8600001 

60.00000 

0.6900000CARD 

E4B 

70.00000 

0.4750000 

80.00000 

0.3400000 

90.00000 

0.2600000CARD 

E4F 

100.0000 

0.2000000 

110.0000 

0.1800000 

120.0000 

0.09000000CARD 

E4G 

130.0000 

0 

.06000000 

140.0000 

0.02500000 

150.0000 

O.OOOOOOCARD 

E4H 

22  DAMPING 

COEFF.  C=35 

;  0  0 

CARD 

El 

0.000000 

1.000000 

0.000000 

0.000000 

1.000000 

CARD 

E2 

0.000000 

35.00000 

0.000000 

0.000000 

0.000000 

O.OOOOOOCARD 

E3 

29  VERY 

STIFF  BELT. 

0  0 

CARD 

El 

0.000000 

-4.000000 

0.000000 

0.000000 

1.000000 

CARD 

E2 

12 

CARD 

E4A 

0.000000 

0.000000 

0.2500000 

4000.000 

0.3333300 

6000.000CARD 

E4B 

0.4166699 

7500.000 

0.5000000 

11640.00 

0.5833300 

14700. OOCARD 

E4C 

0.6666998 

18210.00 

0.7500000 

21600.00 

0.8333300 

25320. OOCARD 

E4D 

0.9166700 

30000.00 

1.000000 

33720.00 

4.000000 

225000. OCARD 

E4E 

31  HARNESS 

N-G  STRAP 

0  0 

CARD 

El 

0.000000 

10.00000 

0.000000 

0.000000 

0.000000 

CARD 

E2 

0.000000 

2500.000 

0.000000 

0.000000 

0.000000 

O.OOOOOOCARD 

E3 

33  BELT 

FRICTION 

0  0 

CARD 

El 

0.000000 

0.000000 

0.9000000 

0.000000 

0.2000000 

CARD 

E2 

34  HAKNESS 

FRICTION 

0  0 

CARD 

El 

0.000000 

0.000000 

1.9900000 

0.000000 

1.9900000 

CARD 

E2 

999 

CARD 

El 

1  RIGHI 

'  SHOULDER  JOINT 

CARD 

E7A 

CARD 

E7B 

-4  12 

CARD 

E7C 

60.07121 

233.6540 

212.3940 

12.03720 

CARD 

E7D 

65.93372 

326.3510 

258.3800 

21.20680 

CARD 

E7D 

82.40853 

356.2500 

217.7140 

11.67850 

CARD 

E7D 

91.21090 

272.0920 

163.5090 

5.066480 

CARD 

E7D 

89.27122 

258.4990 

176.1420 

10.21010 

CARD 

E7D 

89.62075 

288.6960 

176.4290 

6.454900 

CARD 

E7D 

84.04031 

225.0060 

121.3410 

-12.06520 

CARD 

E7D 

80.09302 

195.8620 

95.93260 

-21.24050 

CARD 

E7D 

77.24643 

204.1740 

119.9530 

-11.66900 

CARD 

E7D 

80.60311 

189.2130 

117.3500 

-5.125761 

CARD 

E7D 

99.29332 

188.2180 

75.28892 

-10.24110 

CARD 

E7D 

84.86960 

225.4650 

117.3850 

-6.466969 

CARD 

E7D 

2  LEFT 

SHOULDER  JOINT 

CARD 

E7A 

CARD 

E7B 

-4  12 

CARD 

E7C 

60.07121 

233.6540 

212.3940 

12.03720 

CARD 

E7D 

84.86960 

225.4650 

117.3850 

-6.466969 

CARD 

E7D 

99.29332 

188.2180 

75.28892 

-10.24110 

CARD 

E7D 

80.60311 

189.2130 

117.3500 

-5.125761 

CARD 

E7D 

77.24643 

204.1740 

119.9530 

-11.66900 

CARD 

E7D 

80.09302 

195.8620 

95.93260 

-21,24050 

CARD 

E7D 

84.04031 

225.0060 

121.3410 

-12.06520 

CARD 

E7D 

89.62075 

288.6960 

176.4290 

6.454900 

*  CARD 

E7D 

89.27122 

258.4990 

176.1420 

10.21010 

CARD 

E7D 

91.21090 

272.0920 

163.5090 

5.066480 

CARD 

E7D 

82.40853 

356.2500 

217.7140 

11.67850 

CARD 

E7D 

65.93372 

326.3510 

258.3800 

21.20680 

CARD 

E7D 

4  RIGHT  HIP  JOINT 

CARD 

E7A 

CARD 

E7B 

-4  12 

CARD 

E7C 

63.84941 

826.2141 

676.1950 

-0.2866720 

CARD 

E7D 

58.28361 

826.2141 

676.1950 

-0.2866720 

CARD 

E7D 

38.89900 

826.2141 

676.1950 

-0.2866720 

CARD 

E7D 

36.16820 

826.2141 

676.1950 

-0.2866720 

CARD 

E7D 

39.74170 

826.2141 

676.1950 

-0.2866720 

CARD 

E7D 

50.85709 

826.2141 

676.1950 

-0.2866720 

CARD 

E7D 

63.36500 

826.2141 

676.1950 

-0.2866720 

CARD 

E7D 

47.66820 

826.2141 

676.1950 

-0.2866720 

CARD 

E7D 

37.63620 

826.2141 

676.1950 

-0.2866720 

CARD 

E7D 
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34 

.85600 

826.2141 

676.1950 

38 

.24630 

826.2141 

676.1950 

54 

.14160 

826.2141 

676.1950 

5 

LEFT 

HIP  JOINT 

-4 

12 

63 

.84941 

826.2141 

676.1950 

54 

.14160 

826.2141 

676.1950 

38 

.24630 

826.2141 

676.1950 

34 

.85600 

826.2141 

676.1950 

37 

.63620 

826.2141 

676.1950 

47 

.66820 

826.2141 

676.1950 

63 

.36500 

826.2141 

676.1950 

50 

.85709 

826.2141 

676.1950 

39 

.74170 

826.2141 

676.1950 

36 

.16820 

826.2141 

676.1950 

38 

.89900 

826.2141 

676.1950 

58 

.28361 

826.2141 

676.1950 

8 

RIGHT  ELBOW  JOINT 

-5 

2 

30 

.00000 

90.30700 

271.4830 

34 

.00000 

60.39060 

131.6870 

9 

LEFT 

ELBOW  JOINT 

-5 

2 

30 

.00000 

90.30700 

271.4830 

34 

.00000 

60.39060 

131.6870 

10 

RIGHT  KNEE  JOINT 

-4 

2 

23 

.00000 

52.19941 

441.0020 

34 

.00000 

153.8690 

408.4461 

11 

LEFT 

KNEE  JOINT 

-4 

2 

23 

.00000 

52.19941 

441,0020 

34 

.00000 

153.8690 

408.4461 

12 

RIGHT  ANKLE  JOINT 

-4 

2 

4. 

000000 

179.0880 

67.88540 

4. 

000000 

167.7160 

-12.13370 

15 

LEFT 

ANKLE  JOINT 

-4 

2 

4. 

000000 

179.0880 

67.88540 

4. 

000000 

167.7160 

-12.13370 

999 

3 

5  2 

10  0 

0  0  0 

1 

16  1 

1  13  -20 

-21  0  14 

1 

16  6 

6  13  -20 

-21  0  14 

1 

16  9 

9  13  -20 

-21  0  14 

2 

16  1 

1  13  -20 

-21  0  14 

2 

16  2 

2  13  -20 

-21  0  14 

2 

16  3 

3  13  -20 

-21  0  14 

2 

16  13 

13  13  -22 

-21  0  14 

2 

16  15 

15  13  -22 

-21  0  14 

3 

16  8 

8  13  -22 

-21  0  14 

3 

16  11 

11  13  -22 

-21  0  14 

4 

16  5 

5  13  -22 

-21  0  14 

10 

16  8 

8  13  -22 

-21  0  14 

10 

16  11 

11  13  -22 

-21  0  14 

0 

0  0 

0  0  2 

2  0  1 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

CARD  E7A 

CARD  E7B 

CARD  E7C 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

-0.2866720 

CARD  E7D 

CARD  E7A 

CARD  E7B 

CARD  E7C 

218.8050 

45.72250 

CARD  E7D 

19.23500 

-44.79070 

CARD  E7D 

CARD  E7A 

CARD  E7B 

CARD  E7C 

218.8050 

45.72250 

CARD  E7D 

19.23500 

-44.79070 

CARD  E7D 

CARD  E7A 

CARD  E7B 

CARD  E7C 

176,3340 

CARD  E7D 

250.5860 

CARD  E7D 

CARD  E7A 

CARD  E7B 

CARD  E7C 

176.3340 

CARD  E7D 

250.5860 

CARD  E7D 

CARD  E7A 

CARD  E7B 

CARD  E7C 

133.8270 

CARD  E7D 

132.9940 

CARD  E7D 

CARD  E7A 

CARD  E7B 

CARD  E7C 

133.8270 

CARD  E7D 

132.9940 

CARD  E7D 

CARD  E7A 

2  0  0 

CARD  FIA 

1 

CARD  FIB 

1 

CARD  FIB 

1 

CARD  FIB 

1 

CARD  FIB 

1 

CARD  FIB 

1 

CARD  FIB 

1 

CARD  FIB 

1 

CARD  FIB 

-1 

CARD  FIB 

-1 

CARD  FIB 

1 

CARD  FIB 

1 

CARD  FIB 

1 

CARD  FIB 

110 

10  0 

CARD  F3A 
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6  6  13 

13  3  0 

7  0  19 

CARD  FSB 

6  6  9 

9  3  0 

7  0  19 

CARD  FSB 

7  7  10 

10  3  0 

7  0  19 

CARD  FSB 

7  7  11 

11  3  0 

7  0  19 

CARD  FSB 

9  9  15 

15  3  0 

7  0  19 

CARD  FSB 

10  10  8 

8  3  0 

7  0  19 

CARD  FSB 

11  11  8 

8  3  0 

7  0  19 

CARD  FSB 

13  13  15 

15  3  0 

7  0  19 

CARD  FSB 

0  0  0 

0  0  0 

0  0  0 

0  0  0 

0  0 

CARD  F4A 

4 

CARD  F8A 

9  18  18 

4 

CARD  FSB 

29  0  0 

0  0-0 

.200000  0 

0  0 

CARD  F8C 

16  0  1 

10  0 

0  0  0 

7.000000 

-8.000000 

-8.500000 

F8D1 

0.000000 

0.000000 

0.000000 

2.400000 

22.000000 

-0.300000 

F8D2 

110 

0  0  0 

0  0  33 

2.570000 

-5.650000 

-0.381000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0,000000 

0.000000 

F8D2 

110 

0  0  0 

0  0  33 

3.689000 

-3.594000 

-1.249000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

110 

0  0  0 

0  0  33 

3.812415 

-2.658020 

-1.588378 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0,000000 

0.000000 

F8D2 

110 

0  0  0 

0  0  0 

3.900000 

0.000000 

-2.335000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

1.000000 

0.000000 

F8D2 

110 

0  0  0 

0  0  33 

3.812415 

2.658020 

-1.588378 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

1  10 

0  0  0 

0  0  33 

3.698000 

3.594000 

-1.249000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

110 

0  0  0 

0  0  33 

2.570000 

5.650000 

-0.381000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

16  0  1 

10  0 

0  0  0 

7.000000  ■ 

8.000000 

-8.500000 

F8D1 

0.000000 

0.000000 

0.000000 

2.400000 

22.000000 

-0.300000 

F8D2 

31  0  0 

0  0-0 

.100000  0 

0  0 

CARD  F8C 

16  0  1 

10  0 

0  0  0 

-1.500000 

0.000140 

-32.000000 

F8D1 

0.000000 

0.000000 

0.000000 

0.500000 

0.000000 

0.500000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

-0.592107 

-3.229789 

-6.406097 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

0.197222 

-3.220030 

-6.313524 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

1.000000 

-3.210000 

-6.200000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

1.729085 

-3.200000 

-5.903608 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

2.829251 

-3.125510 

-5.483449 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

3.455790 

-3.042912 

-4.528134 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

3.947432 

-2.836730 

-3.472554 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

4.343934 

-2.664911 

-2.675975 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

4.623846 

-2.493092 

-1.803066 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

4.721993 

-2.286909 

-0.962208 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

4.705504 

-2.149454 

-0.480676 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

4.706552 

-1.875891 

0.730580 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  3  0 

0  0  0 

0  0  34 

4.597802 

-1.756925 

1.387208 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

2  2  0 

0  0  0 

0  0  34 

4.127623 

-0.451147 

-0.532688 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

2  2  0 

0  0  0 

0  0  34 

4.297456 

-0.332015 

0.203587 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

2  2  0 

0  0  0 

0  0  34 

4.098295 

-0.213856 

0.861432 

F8D1 
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1 

31 

16 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

2 

2 

2 

1 

31 

16 

1 

1 

1 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

1  1 

10  0 

0  0  0 

3.900000 

0.000000 

-2.335000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

1.000000 

-2.000000 

F8D2 

0  0 

0  0  ~0 

.100000  0 

0  0 

CARD  F8C 

0  1 

10  0 

0  0  0 

-1.500000 

0.000140 

-32.000000 

F8D1 

0.000000 

0.000000 

0.000000 

0.500000 

0.000000 

0.500000 

F8D2 

3  0 

0  0  0 

0  0  34 

-0.592107 

3.229789 

-6.406097 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

0.197222 

3.220030 

-6.313524 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

1.000000 

3.210000 

-6.200000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

1.729085 

3.200000 

-5.903608 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

2.829251 

3.125510 

-5.483449 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

3.455790 

3.042912 

-4.528134 

F8D1 

0,000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

3.947432 

2.836730 

-3.472554 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

4.343934 

2.664911 

-2.675975 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D 

3  0 

0  0  0 

0  0  34 

4.623846 

2.493092 

-1.803066 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

4.721993 

2.286909 

-0.962208 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

4.705504 

2.149454 

-0.480676 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

4.706552 

1.875891 

0.730580 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

3  0 

0  0  0 

0  0  34 

4.597802 

1.756925 

1.387208 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

2  0 

0  0  0 

0  0  34 

4.127623 

0.451147 

-0.532688 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

2  0 

0  0  0 

0  0  34 

4.297456 

0.332015 

0.203587 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

2  0 

0  0  0 

0  0  34 

4.098295 

0.213856 

0.861432 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

1  1 

10  0 

0  0  0 

3.900000 

0.000000 

-2.335000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

1.000000 

-2.000000 

F8D2 

0  0 

0  0  -0, 

.100000  0 

0  0 

CARD  F8C 

0  1 

10  0 

0  0  0 

18.966637 

-0.000137 

-2.716074 

F8D1 

0.000000 

0.000000 

0.000000 

-1.000000 

0.000000 

1.000000 

F8D2 

1  0 

0  0  0 

0  0  33 

4.478610 

0.000000 

0.868376 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

1  0 

0  0  0 

0  0  33 

4.497387 

-0.000026 

-0.831519 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

F8D2 

1  1 

10  0 

0  0  0 

3.900000 

0.000000 

-2.335000 

F8D1 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

1.000000 

F8D2 

0.0 

0.0 

0.0  0  0 

0  0  0 

CARD  GIA 

13.88821.000010414 
0.000000  13.00000 
0.000000  12.00000 
0.000000  10.00000 
0.000000  5.000000 
0.000000  3.000000 
4.000000  91.24997 
0.000000  53.40000 
0.000000  128.8000 
356.0000  91.25000 
0.000000  53.40000 
0.000000  128.8000 
0.000000  12.00000 
0.000000  70.00000 


-15.12504  0.000000 

0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 


0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 

0.000000  0.000000 


CARD  G2 

321  OCARD  G3A 
32  1  OCARD  G3A 
32  1  OCARD  G3A 
321  OCARD  G3A 
321  OCARD  G3A 
321  OCARD  G3A 
321  OCARD  G3A 
321  OCARD  G3A 
321  OCARD  G3A 
321  OCARD  G3A 
321  OCARD  G3A 
321  OCARD  G3A 
321  12CARD  G3A 
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Portion  of  Example.aou  File 


Example.aou  is  the  primary  output  file  (Unit  6)  of  the  sled  test  simulation.  Because  of  its  large  size, 
only  a  portion  of  the  file  (results  up  to  10  msec)  is  presented  here.  The  first  part  of  the  file  contains 
the  input  data  with  all  the  variables  clearly  labeled.  The  second  part  of  the  file  contains  the  simulation 
results  at  predetermined  time  intervals. 


58 


AAMRL  ARTICULATED  TOTAL  BODY  (ATB)  MODEL  PAGE 
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SIMULATION  OF  THE  HUMAN  VOLUNTEER  SLED  TEST 
USING  NEWLY  DEVELOPED  HUMAN  JOINT  CHARACTERISTICS 
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FIRST  PART  OF  FUNCTION  -  8  TABULAR  POINTS 
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FUNCTION  NO.  29  VERY  STIFF  BELT.  NTI(29)  =  126  CARDS  E 


FIRST  PART  OF  FUNCTION  -  12  TABULAR  POINTS 
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FUNCTION  NO.  34  HARNESS  FRICTION  NTI(34)  =  172  CARDS 


FUNCTION  IS  CONSTANT  1.990000 

1  JOINT  FORCE  FUNCTION  NO.  1  RIGHT  SHOULDER  JOINT  NTI  (  1)  =  177  PAGE 

CARDS  E.7 
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-180.00  23.000  52.19941  441.0020  176.3340 

0.00  34.000  153.8690  408.4461  250.5860 

1  JOINT  FORCE  FUNCTION  NO.  12  RIGHT  ANKLE  JOINT  NTI(12)  =461  PAGE 
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Example  Time  History  Files 


Time  history  files  are  those  files  with  output  logical  units  greater  than  21 .  The  data  in  time  history 
files  are  arranged  in  columns.  These  column  data  can  be  easily  ported  to  any  spreadsheet  or  graphics 
software  for  further  result  analysis.  There  are  a  total  of  33  time  history  files,  from  example.t21  up 
to  example.t53,  being  generated  by  using  example.ain  as  the  ATB  input  file.  In  this  section,  only 
one  file  for  each  type  of  time  history  output  is  presented.  The  included  files  are: 

1 .  Example.t21 :  Point  linear  accelerations. 

2.  Example.t22:  Point  linear  velocities. 

3.  Example.t23:  Point  linear  positions. 

4.  Example.t24:  Point  angular  accelerations. 

5.  Example.t25:  Point  angular  velocities. 

6.  Example.t26:  Point  angular  positions. 

7.  Example.t27:  Joint  parameters. 

8.  Example.t37:  Joint  forces  and  torques. 

9.  Example.t38:  Plane/segment  contacts. 

10.  Example.t45;  Harness  belt  contacts. 

11.  Example.t47:  Ellipsoid/ellipsoid  contacts. 

Each  file’ s  output  unit  is  shown  by  the  highlighted  file  page  number. 
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APPENDIX  B 


Developing  New  Deformable  Segment  Models 

There  are  three  steps  to  creating  a  deformable  body  model  for  the  ATB: 

1)  Develop  a  finite  element  model  of  the  segment  and  perform  modal  analysis  to  determine  the  first 
two  natural  frequencies  and  their  corresponding  mode  shapes. 

2)  Write  an  ASCII  file  containing  the  nodal  positions,  nodal  masses,  frequencies,  and  mode  shapes. 
The  user  must  make  sure  that  this  file  has  the  format  described  in  section  B.2. 

3)  Modify  the  ATB  input  file  according  to  the  ATB  Input  Description. 

B.  1  Segment  Finite  Element  Model 

A  finite  element  model  of  the  segment  must  be  created  and  modal  analysis  must  be  performed  to 
determine  its  natural  frequencies  and  mode  shapes.  The  required  information  are  number  modes, 
number  of  nodes,  nodal  positions,  nodal  masses,  frequencies,  and  mode  shapes. 

B.2  Deformable  Body  Input  Data  File 

The  finite  element  modal  analysis  information  must  be  written  in  an  unformatted  ASCII  file  in  the 
following  order: 

1)  number  of  nodes  (NNOD)  and  number  of  modes  (NMOD);  No.  of  data  =  2 

2)  nodal  positions  with  respect  to  the  body  reference  frame;  No.  of  data  =  3*NNOD 

3)  nodal  masses;  No.  of  data  =  NNOD 

4)  natural  frequencies  in  Hz;  No.  of  data  =  NMOD 

5)  mode  shapes  (eigenvectors);  No.  of  data=  6*NNOD*NMOD 

This  procedure  has  been  simplified  when  ANSYS®  or  ALGOR®  is  used.  A  FORTRAN  program 
named  “atbalgor.for”  has  been  written  which  creates  the  required  ASCII  file  and  assigns  the  name 
“filename.dat”  using  the  “filename.l”  &  “filename.frq”  files  produced  by  ALGOR.  Another 
FORTRAN  program  named  “atbansys.fof  ’  has  also  been  written  for  ANSYS  which  also  creates  the 
same  ASCII  file.  These  programs  are  included  with  the  ATB  model. 
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